Avoided Phase Transitions and Glassy Dynamics In Geometrically Frustrated 

Systems and Non-Abelian Theories 
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We demonstrate that the application of any external uniform non-abelian gauge background, no 
matter how small, leads to a greatly enhanced degeneracy. This degeneracy is so large that even a 
non-abelian background field of infinitesimal strength leads to a shocking change in the thermody- 
namics. The critical temperature might be discontinuously depressed and an "avoided critical point" 
will emerge. We focus on how this arises in models previously employed to describe the microscopies 
of metallic glasses and correctly predicted the structure factor peaks. Some of the best fits, to date, 
to the dynamics of supercooled liquids were inspired by such notions for which we now provide 
a suggestive microscopic basis. We generalize the Mermin- Wagner inequality to high dimensions 
and discuss how extensive configurational entropy may be computed, by replica calculations, for a 
multitude of glass models (including non-Abelian gauge backgrounds). This extensive configura- 
tional entropy then allows a possible derivation of Vogel-Fulcher dynamics. We fortify earlier ideas 
suggesting avoided critical dynamics. 



I. INTRODUCTION 

Geometric frustration [1] spans a wide range of 
systems- Frank-Kasper phases [2] [3] , [4] , disordered ma- 
terials, clathrates, and complex fluids such as the chol- 
steric blue phases [5], and amphiphilic membranes [6]. 
Perhaps most notable of these materials are metallic 
glasses, and possibly glasses in general, whose ideal icosa- 
hedral (and more complicated extensions) packing is not 
compatible with EucUdean space filling. Such packings 
have seen a revival in the study of quasicrystals. Till now 
many theoretical treatments, e.g. [7], have treated geo- 
metrical frustration as a perturbative effect akin to the 
smooth suppression of the superconductivity by the in- 
sertion of a frustrating magnetic field and the appearance 
of topological defects (vortices) in type II superconduc- 
tors. A first hint that the physics might be far richer and 
reward us with interesting surprises was found in [10], [11] 
where the insertion of an infinite range Coulomb inter- 
action was shown to be non-perturbative. By analogy to 
the 1 /r gravitational potentials generated by geometrical 
deformations, all of this suggested that a similar occur- 
rence might arise in geometrically frustrated systems. 

Here we will show how an explicit infinitesimal geomet- 
ric frustration may lead to a dramatic, non-perturbative, 
effect on the thermodynamics in a multitude of systems. 
As a result of geometric frustration, the critical tempera- 
ture might be discontinuously depressed and an "avoided 
critical point" will emerge. The physical origin of the 
discontinuity is a greatly enhanced degeneracy: the geo- 
metric frustration thwarts ideal order (or crystallization) 
and leads to many competing low energy states. In this 
article we will focus on how this arises in various 0{n) 
systems immersed in uniform SO(n) backgrounds. Sim- 
ilar models [5], [7] were previously employed to describe 
the microscopies of metallic glasses and correctly pre- 



dicted the locations of the structure factor peaks [7] . We 
will show that that such models exhibit "avoided critical 
behavior" : a genuine true non-analyticity only at very 
low temperatures as compared to the (avoided) critical 
temperature of the system sans any frustrating external 
non-abehan gauge background. Some of the best fits [12], 
to date, to the dynamics of supercooled liquids were in- 
spired by such notions for which we now provide a firmer 
microscopic basis. We will compare one of the fitted pa- 
rameters in these plots (the magnitude of the "avoided 
critical temperature" for each of the supercooled liquids) 
and demonstrate that, on average, these might be well 
predicted by our microscopic picture. 

To couch geometrical frustration in a general field theo- 
retical setting, we will consider systems whose continuum 
Lagrangian density may be written in the form 
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with a covariant derivative D^[x) — — i6A^[x) appro- 
priate for parallel transport on a surface whose curvature 
is (l/R) = d. As we will detail below, the strength of the 
coupling to the background gauge field [0) sets the degree 
of geometric frustration present in the system. 

Let us quickly sketch what unfolds as the frustration 
is introduced. Matters become most transparent in mo- 
mentum space which diagonalizes the quadratic part of 
the action. When — the quadratic part of the ac- 
tion has a single minimum at wavenumber k = Q. In this 
standard case, long wavelength (small |fc|) fluctuations 
can destroy order in dimensions d < 2 and indeed they 
do so for continuous flelds (f). In high dimensions {d > 2), 
order reasserts itself at temperatures T < Tc — 0{1). 

This picture undergoes dramatic ramifications as the 
frustrating field is switched on. As seen, in Fig.(l), the 
minimum in k space is no longer point like but rather 



1 



extends over a (d-1) dimensional manifold whose radius 
is set by the geometrical frustration 9. For any frustra- 
tion 6, no matter how small, the manifold of minimizing 
modes M is a high dimensional object spanning many 
modes. All Fourier modes on this surface are degenerate 
leading to an exceptionally high configurational degen- 
eracy in real space. In the aftermath, an effective one 
dimensional behavior results (which can be made precise 
in the n — > oo limit) with no ordering transition in sight. 
Only at the unfrustrated point 9 = docs the minimiz- 
ing manifold M shrink to the origin, and our old intuition 
regarding the long wavelength fluctuations is regained. 

For simplicity, we will mostly consider a translationally 
invariant with non-commuting uniform gauge connec- 
tions Ajj,. Our results can be extended to systems with 
arbitrary non-uniform D^{x) in which case our analysis 
may be repeated with the Fourier basis index replaced 
by an index labeling the more complicated diagonaliz- 
ing eigenbasis (the minimizing manifold M will merely 
appear as an identical object in this basis). 

We will consider non-Abelian gauge background fields 
{Afj_{x)} having no dynamics of their own and merely en- 
forcing geometrical frustration. To make the gauge field 
a true dynamical variable, we will need to add terms 
involving derivatives; the simplest standard gauge in- 
variant choice is the squared field (-F)ti/)^ where the field 
Fij,v = [Dix,D^]. The resulting well known Lagrangian 

r _ i P PC*^ 

■^matter -I- gauge ~ 4 

+l\D,r\'-^\$\' + ^_\$\' + .: (2) 

Theories involving only the first term are frequently re- 
ferred to as "pure gauge theories" . We will prove that in 
the opposite extreme of Eqn.(l) (non-dynamic or small 
j4ju) a new surprising "avoided critical behavior" may oc- 
cur: the effect of coupling to the non-Abelian gauges 
is non-perturbative. The theory suddenly becomes ex- 
tremely sensitive to thermal effects. The density of states 
becomes, to a certain extent, one-dimensional. This also 
applies to some more general cases in which the effective 
action after integrating out the gauge fields will be of a 
form similar to Eqn.(l). For instance, if the kinetic term 
in the gauge fields is replaced by {F^^ — f^v)"^ with speci- 
fied frustrations /^^ then after integrating out the gauge 
fields An{x) the resulting action density in the matter 
fields (f) will exhibit frustration. Such a theoretical pre- 
diction has serious implications for many geometrically 
frustrated systems which will form a new universality 
class of their own. As we will aim to highlight in the pa- 
per, the very unusual "glass transition" might be linked 
at its very core to the non-perturbative physics spawned 
by non-Abelian background fields (although of a more 
complicated form than of a fixed uniform background). 

Our main conclusion in the continuum limit is 
schematically summarized in the large n phase diagram 




FIG. 1. A schematic of the interaction kernel, v{k), minima 
in Fourier space. In all of the systems that we consider here, 
the relevant continuum limit v{k) attains its minima on shells 
\k\ = q > Q. In the presence of lattice point group symme- 
try terms, the rotational symmetry is lifted and the minima 
appear at a discrete number of points (labeled by "X" in the 
figure) . 

of Fig. (2). Here, the horizontal axis denotes the strength 
6 = of a uniform non-abelian SO{n) gauge back- 
ground. The vertical axis denotes the critical temper- 
ature of a three dimensional 0{n) spin system. In the 
absence of a uniform non-Abelian background field (i.e. 
at 9 = 0), the spin system undergoes a phase transition 
at a finite T = T'c(O) which is order of the exchange con- 
stant J. When a non-abelian background field is applied, 
the system is unable to order at finite temperatures. 

We also report on a generalized Mermin- Wagner in- 
equality which allows us to connect the magnitude of the 
order parameter (the absence of cntropic effects) to the 
characteristic relaxation times present in any system (in- 
cluding glass models). We will also show how to derive 
Vogel-Fulchcr dynamics for a multitude of glass models 
(including non-Abelian gauge backgrounds), and try to 
fortiiy earlier ideas suggesting avoided critical dynamics. 

II. DEFINITIONS AND CONVENTIONS 

Throughout, we will examine the interaction kernel v 
appearing in the quadratic part of the Hamiltonian 

^=\j -^^z'^^i^^)^^^^)^^^-^)- (3) 

This is the Fourier transform of the real space H = 
^ / dfix J df'yVij{x — y)(j)i{x)(j)j{y). For spins lattice sys- 
tems we will consider Hamiltonians of the type 

H=\Yy{x-y)[S{x)-Sm- (4) 

The sites x and y lie on a (generally hypercubic) lat- 
tice of size iV, and {Si{x)} denote the i-th component 
of the spin S situated at a certain site x. The normal- 
ized spins have n components. Throughout, we employ 
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FIG. 2. The typical phase diagram encountered in n > 2 

realizations of our models. When the system is unfrustrated, 
the critical temperature in dimensions d > 2 is finite (and of 
order 0(1)). When the minimizing manifold M is a (d — 1) 
dimensional shell, the critical temperature drops to zero. In 
the above, 1 /R corresponds to the inverse radius of the sphere 
on which ideal order (icosahcdral or other) can prevail. 1/R 
scales as q- the radius of the minimizing shell M is A; space. 



the non-symmetrical Fourier basis convention (/(fc) = 
EsF(f)e-^^-; F{x) = j^Y.kfik^'^-^) ■ The vec- 
tor q denotes a wave- vector minimizing the kernel v{k): 
v{q) = mink{v{k)}. In the usual ferromagnetic models 
^ = is the only minimizing mode. For simplicity, we set 
the lattice constant to unity and consider a lattice (of size 
L) with periodic boundary conditions. The wave-vector 
components ki = where ri is an integer (and the real 
space coordinates .t; are integers) . In all systems that we 
study in this paper the Fourier modes gi minimizing the 
interaction kernel v{k) lie on a (d-1) dimensional shell. 
We often refer to this high dimensional minimizing man- 
ifold as M. The high dimensional k space redundancy 
leads to a complicated "energy landscape" for real space 
configurations. As there are many zero modes tangent to 
M that cost no energy, the system is very fragile to per- 
turbations and in the large n limit. To = 0. On a lattice, 
the continuous rotational symmetry is lifted. Here, we 
find a finite number of minimizing modes (denoted by 
an "X" in Fig.(l)). Notwithstanding, the critical tem- 
perature will typically be much lower than that of the 
unfrustrated system with a single minimum at fc = 0. 
The critical temperature will drop discontinuously (from 
Tc = (9(1) to Tc «; 0(1)) as M expands from the origin. 



III. ICOSAHEDRAL ORDER IN BRIEF 

Unlike other sections to follow, all presented here is 
a review of earlier work [1]. [5], [7] which inspired the 
current publication. When a liquid is supercooled (when 
cooled so rapidly that it may not crystallize), and cooled 
further down until its relaxation times become effectively 
infinite it becomes a "glass" . In slow solidification into an 
ordered crystal, the atoms/molecules readjust their posi- 



tions, slowly "computing" their ideal crystalline positions 
and veer toward them. Now, consider the situation, for 
the supercooled liquid, when the "computation" needs to 
be done almost on the spot. Here, the atoms/molecules 
can only probe their immediate surrounding to minimize 
energy costs of bonds with their nearest neighbors. For 
a Lennard- Jones interaction, each atom would like to be 
some distance rmin (the distance at which the Lennard 
Jones like potential V{r) is minimized) away from all of 
its neighbors in order to minimize each bond individu- 
ally. Here, crystallization is trivial in two dimensions- 
the local minimum can be extended to tile the entire 
plane. Locally and globally, the ideal packing is hexago- 
nal with each atom surrounded by its six nearest neigh- 
bors all equidistant from it at the ideal distance Tmin- 
The local minimum that the atoms may find is also the 
global minimum. Such a fortunate occurrence does not 
arise in three dimensions. Four atoms may be equidis- 
tant from each other to form the vertices of a perfect 
tetrahedron. However, around a given edge of the tetra- 
hedron we may not pack five other tetrahedrons perfectly: 
We can almost do so leaving a 7° void. Similarly, parti- 
cles whose locally preferred structure is an icosahedron 
(of 13 particles) cannot globally tile all of space as seen 
by the "illegal" icosahcdral five-fold rotational symmetry 
[16]. If this is correct, then microscopically, local clusters 
should appear. Certain supercooled metallic glasses of- 
fer a good testing ground for these notions. There is 
evidence that short-range order in under-cooled liquids 
and metallic glasses is icosahcdral [8] . A theoretical con- 
struct invented long ago [1,3,5] is a high dimensional ref- 
erence system in which perfect local icosahcdral order, 
with equal bond distances, could easily tile large regions 
without any voids. Endowing three dimensional space 
with a small curvature, the 7° void can be compressed 
to zero: Icosahcdral clustering is perfect in a slightly 
curved three dimensional space. An ideal, icosahcdral 
crystal on the surface of a 4-sphcre (polytope "{3, 3, 5}" ), 
can be found. It consists of 120 particles embedded on 
the slightly curved spherical surface (the three dimen- 
sional boundary of a sphere in four dimensions). We ( 
[12]) have taken this notion to the extreme and argued 
that the dynamics of glasses is controlled by the vicinity 
to the melting transition of ideal crystal that may form 
on the curved sphere. Regions of short-range {3,3,5} 
order in are broken up by an array of —72° disclination 
lines, forced in by "frustration"- the incompatibility of 
flat space with a space filling icosahcdral crystal. The 
Frank-Kasper phases of transition metal alloys are or- 
dered arrays of disclination lines in such an icosahcdral 
medium. The order parameters <5„ ,„a, mi, arc obtained by 
projecting the local particle configuration onto the hyper- 
spherical harmonics l^,ma,m6 (the spherical harmonics 
("yj'""-s) for a 4-sphere) of a tangent four-sphere which 
can accommodate icosahcdral order [7] . 

In three dimensions, a 3-sphere is attached to every 



3 




X-axis 



^AVERAGING ^ 
VOLUME 

FIG. 3. A one dimensional analogue. We project the mass 

density from the real-space .t— axis onto a disk of radius R 
(with internal coordinate 9) associated with each point on the 
X-axis and expand in harmonics {p{x,9) = pm(x)e""''). 
Ideally, the mass density on the line matches that ob- 
tained by rolling the disk and the local energy functional 
[dx — de/R]p{:i:,9)]'^ is minimized. Integrating over 9 we find 
the one dimensional analogue of the 50(4) action below. 

point X in real space. The mass density is projected onto 

these spheres. The density at all x on the tangent sphere 
(u is the internal coordinate on the sphere), 

P{X,U)= Qn;m+,m-{x)Yn;m+,m-{u), (5) 

n;m+ ,m_ 

Upon expansion in the hyper-spherical harmonics 

^^;m+,m_* The 771^ SUbscriptS in ^^;m-|-,m_ 

refer to the 

two SU{2) representations (of generators A± of spin n/2) 
that the 5*0(4) representation can be decomposed into. 
The Hamiltonian 

H =lj2i^-\(^--'^^oa)Qn\'+f^n\Q\% (6) 
n 

where we organize the {Qn-.m+m-} multipoles for 
each n into (n + 1)^ dimensional vectors in the 
|n/2, m_|_; ri/2, m_) basis. In a large n analysis, the sec- 
ond (mass) term portrays tlic effect of normalization. 
Fourier transforming and diagonalizing [7], 

^=1 H |a^,,(fc)|^[/^„A^,,(fc)+Mn], (7) 

k,M,i,n 

with M the eigenvalues of k- [A++A^], the amplitudes of 

the i-th eigcnmodc af(fc), and {A"| j(fc)} the eigenmode 
spectrum. Diagonalization yields, in all cases, a set of 
eigenvalues {A^(|fc|)} with the global minimum (for each 
n > 0) occurring at = (7^ > 0, in agreement with struc- 
ture factor measurements (displaying a series of peaks at 
these minimizing wave- vectors). By rotational symme- 
try, the eigenvalues depend on the magnitude \k\ of the 
wavenumber. Thus, as the minima do not occur at the 
origin, shells of minimizing modes (as in Fig.(l)) occur. 
We now define 

min{i^„Aj.(|^|)} = i^„„.„Ajrj„(g), (8) 



where rimin and ?'ijii„ denote the values of these indices at 
which the global minimum is attained- the correspond- 
ing minimizing wavenumber is denoted by q. By ro- 
tational symmetry, all that matters is the modulus q. 
As expected, q scales with the geometric frustration 9. 
The values of n allowed by icosahedral symmetry are 
n = 0, 12, 20, 24, 30, 32, 36. The smallest non-trivial n 
corresponds to a (n-|- 1)'^ = 169 complex component vec- 
tor Qn,Tn-\-,7n- 

_ . For a discussion on how these ideas may 
also be applied to wave functions on a sphere [14]. 

IV. NON-PERTURBATIVE THERMODYNAMICS 

We will argue, via a large n analysis, a Mermin- Wagner 
inequality for similar models, and a fluctuation analysis 
[24], [33], that these and other continuum 

glass models do not order at finite temperatures. 

Any finite frustration 6, no matter how small, drives 
these systems away from criticality at finite tempera- 
tures. In many solvable cases, this non-perturbative ef- 
fect has nothing to do with the breakdown of pertur- 
bation theory due to long range interactions. For any 
value of the frustration 0, the minimizing manifold of 
modes in Fourier space is high dimensional. The high 
dimensionality allows for the proliferation of many zero 
energy modes which easily erase any sign of potential or- 
der. When small point group symmetry terms lift the 
rotational symmetry of continuum models, the critical 
temperature becomes finite yet still minute as compared 
to the ordering temperature of the unfrustrated system, 

o{i) = Tc{0 = 0) » me + 0). (9) 



V. THE LARGE N PHASE DIAGRAM 

Whenever the minimizing mode q is finite and a (d- 
1) shell of minimizing modes occurs (as in Fig.(l)), a 
large n calculation reveals that the critical temperature, 
rps-pherical ^ q. As a function of the radius of the min- 
imizing shell M in Fig.(l), this drop occurs discontinu- 
ously. 

In the large n (or spherical) limit, we subject the sys- 
tem to a normalization constraint of the (n -|- 1)^ (that 
of n = nmin) component vector Q, 

E I«mT(^)I' = i- (10) 

M,i,k 

We solve the system of Eqn.(7), subject to the latter 
(spherical or large n) constraint, to find that the critical 
temperature. 
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<pk 



(27r)3 



(11) 



As all terms in the integrand are explicitly positive, we 
may bound the right-hand side by its contribution from 
(imin, M = 0). If the (imin, M = 0) eigenvalue is analytic 
near its minimum at |fc| = then using the shorthand 

i^(|fc|) = -fi^ramin-^o'imi„(^)' ^ trivial bound (from the re- 
gion [q — 5,q + S\) reads 



> 



{q-5f 



dk 



q-5 



[F"{q)]{k-q)^ 



(12) 



The integral in Eqn.(12) diverges. Here we assumed that 
the second derivative of F is finite at A: = g; if the second 
derivative vanishes then the divergence of the right-hand 
side of Eqn.(ll) is even more dramatic. The divergence 
signals that the critical temperature Tc ~ 0. Critics might 
argue that this might be an artifact of the large n approx- 
imation and that, "for well known" reasons, this or that 
may happen (e.g. fluctuation driven first order phase 
transitions [35]). To counter such arguments, we prove 
in a later section a generalized Mermin- Wagner theorem 
of related rotationally symmetric spin models disallowing 
magnon dispersion relations about any ordered system at 
arbitrarily low temperatures T = 0+! Even if a fluctua- 
tion driven first order transition occurs, as we show (by 
extending replica calculations), the system may be frozen 
in a glass state before reaching the transition, so its exis- 
tence might be immaterial. For the moment, we simply 
note that the lowest lying vector (5n=i2 bas 169 complex 
components. This number is large and a 1/n expansion 
seems plausible. In contrast to certain expectations, it 
is seen that many beautiful early constructs [1,7,5] pre- 
dict a vanishing Tc for any finite non-Abelian coupling if 
no symmetry breaking terms arc present. As shown in 
Fig. (2), the phase diagram corresponding to these ideal 
metallic glasses, and plausibly to dense random packed 
supercooled liquids in general, shows an "avoided critical 
point" - a spike for tunable g = in the interaction kernel, 
a genuine non-analyticity in the (thermo)dynamic func- 
tions at finite coupling to the uniform gauge background 
only at = 0. 





FIG. 4. The self energy corrections. The thin dashed lines 
denote bare interactions. Thick dashed lines represent dressed 
interactions (i.e. the bare interactions screened by a geometric 
series of bubble diagrams). The solid lines denote propaga- 
tors. E° is the 0(1) parquet diagram. E^^^'^^' are 0(l/n) 
corrections. 



Go\k)=Kr,„ 



\nmi 



(k)+ro. 



(14) 



The zeroth order self-energy is a momentum independent 

constant 



-J 



_fk_ 

(27r)3 



Gi. 



Sk). 



(15) 



The zeroth order self-energy contributions from all 
other {non-{nmin,imin)) eigenvalues are strictly posi- 
tive. When the mass gap vanishes, i.e. r^nin = 
—Kn^in^o'i^.^iq), the integral in Eqn.(15) diverges. By 
the Dyson equation, 



G-\k)=Go\k) + i:{k). 



(16) 



A divergent 5]^°) ~ (r — rmin)~^^^ signals the inabil- 
ity of the mass gap to vanish. Only at ro = — oo can 
r = — -f^'rinjin^oT'" (*?)• the conventional (finite Tc) 
case, the bare mass rg = a(T — Tq) codes for a linear 
temperature scale. The divergence that we obtain here 
implies that = and that consequently a linear tem- 
perature scale is void in its environs. A similar analysis 
for the quantum case is detailed in the appendix. There 
arc two first order 1/n diagrams and, after regrouping 
diagrams self-consistently, six important second order di- 
agrams). By re-summing [24] diagrams self-consistently, 
the zeroth order divergence forbidding finite temperature 
ordering may stay in tact. 



VII. NEW SO{N) BACKGROUND SPIN MODELS 



VI. AVOIDED CRITIC ALITY VIA 1/N 

Within a 1/n diagrammatic framework, the absence of 
a finite temperature phase transition, in the large n limit, 
is seen by the divergence of the lowest order (tadpole) 
contribution to the self energy (Fig. 4). With the self- 
consistent propagators, 

G7l,n^Jk) = K^_\---^{k) + r (13) 

for the lowest energy mode {nmimimin) of relevance at 
low temperatures, with the bare propagator 



Here we, by inventing a simple "soccer-ball model", we 
intuitively illustrate a simple basic premise: 

A non-Abelian background leads to High Dimensional 
minimizing shells in k-space. 

The radius of the shell is set by the strength of the 
non-Abelian background. The profound increase in the 
dimensionality of the minimizing modes leads to a highly 
non-perturbative effect on the thermodynamics. In mo- 
mentum space, there are many zero-energy, "Goldstone" 
like, modes tangent to the minimizing manifold. These 
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fluctuations lead to a vanishing spin stiffness. All of this 
is in stark contrast to the innocuous effect of abelian 
gauge backgrounds: A U{1) minimal substitution v{k) = 
{k — eA)"^, with a constant A, trivially shifts the point 
minimum in - space. For non- Abelian coupling the ef- 
fect is much more subtle: the minimizing m,anifold is a 
shell spanning an infinite number of Fourier points. As 
the original models of [7] only allow for representations of 
SO (3) or SO (4) with a high number of components (large 
n indeed) for which calculational results may not be too 
transparent, we will now invent a simple toy spin model 
to gain intuition: the "soccer-ball model" . The interplay 
between the external spatial and internal spin direction 
that it possesses is present in very few other spin systems 
of a very different character (e.g. orbital Jahn- Teller and 
Kugcl-Khomskii like compass models). Consider three 
component (Heisenberg) spins on the square lattice with 
SO (3) couplings: 



H = -Y,S{x)*S{y). 



(17) 



{x,y) 



Here the * operation denotes to the scalar product be- 
tween a spin S{x) attached to the south pole a sphere (a 
"soccer-ball") and subsequently rolled to the neighbor- 
ing site y and the spin S{y) situated at y ; the action is 
locally minimized by a spin rotation in the [f.n] plane for 
a lattice displacement along the £ — axis. For a rolling 
angle 9, we will now consider a spin configuration having 
as its only non-vanishing Fourier components S'(±^ei), 
namely, the configuration 



(S'i,S'2) = (cosfei, sintei) 



(18) 



with S'l and 5*2 the x and y components of the spin and 
Xi the spatial x coordinate. This spiral configuration will 
saturate all bonds along e\, and satisfy S{x) = S{x -\- 62) 
(just as for the uniform (fc = 0) configuration) . The en- 
ergy of the above configuration is lower than that of a 
fc = mode. In the continuum limit, the kernel is in- 
variant under rotations of k. As the minimum does not 
occur at the origin, a shell of minimizing modes is seen 
to exist (a {d — w) dimensional manifold of minimizing 
modes is generated for 0{n = d + 2 — w) spins with the 
above "rolling" SO{n) action for all but {w — 1) direc- 
tions assigned ferromagnetic couplings). Alternatively, 
in the continuum limit, the differential action evaluated 
for the above configuration along any ray not parallel to 
62 is lower than that for the uniform configuration. The 
continuum Hamiltonian 



H = 



1 



ID^S^l^d^x 



S^{k)vij{k)Sji-k)d'^k, 



(19) 



a nonlinear a model on a curved surface. The covariant 
derivative 



(20) 



where is the v component of the ^ = 1 angular mo- 
mentum generator (ei2 = —£21 = 1)- Demanding that 
the "rolled" vector increment be real suggests the an- 
gular momentum representation (Lj)j7; = —teijk- 

The 

interaction kernel pursues the form: 

v,^j{k) = k'^S.^j + 2ee^,k^{L,),^, + e^{Ll + Ll),^j. (21) 

The eigenvalues are invariant under planar rotations of 
k. The minima occur on the ((2-1) dimensional) shell: 



For the lattice version Vij{k) 



(22) 
±''=^i?^(±6») 



where Re{±d) is the rolling matrix (rotation matrix 
about the t ^ ^-axis ). The continuum representation 
Li = —ieijk (for the change of the spin 6S = 6(f) x S 
incurred by a rotation of angle S(j)) emerges trivially. 

Exploiting rotational invariance in the continuum 
limit, we may always rotate a pure mode q of unit (0(n)) 
spins such that it lies parallel to ei. All bonds parallel 
to 62 will entail the same energy penalty as the k = 
mode and to minimize energy cost along ei we simply 
set q = 9ei. We now relax the constraints of locally 
normalized real spins and allow the spins to be complex 
and only satisfy a global normalization constraint. Any 
ground states found here will obviously be of equal or 
lower energy than those of three-component Heisenberg 
spins. In a spherical model of complex spins with global 
normalization: '^g\S{x)\'^ = N enforced by a Lagrange 
multiplier /i, the continuum limit kernel reads 

+ +11 2l0ky 

k'^ + e'^ + fi 2iek^ 

-2i6ky -2idk^ k'^ +29"^ + fj,^ 

Due to an inherent coupling between real and imaginary 
components, |g| = ^ may not be enforced. 

As they must, all eigenvalues degenerate onto the 
canonical fc^ dispersion in the limit 0—^0. Two of the 
eigenvalues 



A'± = A± 



x'^ = Xo + ^I = k^ + 

■k^ + -0^ ±-^/l6k^02 

2 2 



(23) 



min{A'_} 

k 



(24) 



are non-analytic in The global minimum 

occurs on the (2-1) dimensional annulus |g| = ^^^|^|- 
The corresponding eigenvector is 

-2i0kx 



\u-) = AA|1,. 



-2i0L, 



-6*2/2 



i6iVl6/c2 + 6i2' 



(25) 
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where the R,I superscripts denote real/imaginary com- 
ponents respectively. The minimum of the non- 
oscillatory "ferromagnetic" mode occurs at the origin , 

min{A;} = X'^{k = 0) = ij + e\ (26) 

k 

A Landau-Ginzburg type analysis ( [24], [33]) shows 
strong (non-logarithmic) divergent fluctuations from k 
q in directions (within the 6-dimensional parameter 
space (sf s^'^, s^'^) ) transverse to the eigenvector of 
Eqn.(25) (such that the norm |s| does not change to low- 
est order) . The generalization of our "soccer-ball" model 
to SO (4) couplings on a three dimensional lattice is dis- 
cussed in Appendix(B). The reader might be dismayed 
that here we analyzed a situation in which the spins are 
complex. Here we attempt to capture as best as possi- 
ble the complicated non-Abelian theories of glasses. In 
the actual high dimensional multipole expansion as in 
Eqs.(6) the "spin components" Q are complex. 

VIII. OTHER INTERACTIONS 

Non-Abelian background fields are not the only poten- 
tial sources of glass formation and amorphous "clump" 
configurations. Many more physical Euclidean space 
potentials seem to do the trick. All of these interac- 
tions seem to share the generic feature that in Fourier 
space- the minimizing modes form a high dimensional 
shell (Fig.(l)) leading to an extremely high degeneracy 
in real space. This in turn leads to a "complex energy 
landscape" with many competing minima. 

A. The Coulomb Frustrated Ferromagnet 

My collaborators and I have [12] modeled the geomet- 
ric frustration felt by molecules as they attempt to pack 
into a locally preferred structure by a long-range coulomb 
interaction forbidding the ideal local structure to tile all 
of space. In the continuum limit, the Hamiltonian 

+ j (fx j d^x'(l){x)V{x- x')(l){x'), (27) 

with the potential V{r) = e^/r. The minimizing modes 

occur on shells as in Fig.(l) with \q\ ~ e^/^. Here we will 
derive such a model from the underlying microscopies. 

B. Kac-type Potentials 

These two-body interactions, 



V{\x\) = yo7''<^(7l^l) (28) 

with the step function (j){y) = 1 for y < 1 and zero 
otherwise and ¥0 = 0^ [22]. The constant a controls 
the integral strength of the potential and 7 dictates the 
range of the potential. We note that, physically, this 
potential corresponds to a weak core interaction inside 
a sphere of radius 7"^. These interactions may model 
glass formation [19]. The Fourier transformed interac- 
tion is minimized at finite wave-numbers (as easily seen, 
it has a maximum at A; = when all the Fourier phases 
add coherently). The Fourier transform of Eqn.(28), 

v{\k\) = 47rVb7-''[sin(fc/7) - (k/-/) cos{k /-/)]/ k^. (29) 

This kernel displays its minimum at = q > with 
q/j 5.76. The minimum indeed occurs on a (d—l) di- 
mensional shell such as that drawn in Fig.(l). The kernel 
v{k) is analytic near |fc| = g allowing us to Taylor expand 
about the minimum. By superposing Kac-type poten- 
tials, we may generate discrete, "digitized", versions of 
the Lennard- Jones potential. For instance, 

V{x)=J2VnUhi\^l (30) 

with Vb = » —Vi > 0, 71 = 1/ri, and 72 = l/r2 
corresponds to a potential V{x) which is (Vq + Vi)^0 
inside a sphere of radius ri and is Vi < for ri < r < r2, 
and vanishes for all r > r2- Superposing the Fourier 
transformed kernels of Eqn.(29), shows a minimum on 
a spherical shell of finite radius \k\ = q > 0. Although 
the absolute minimum of a particles interacting via a 
Lennard-Jones is an FCC crystal, whenever the liquid is 
cooled rapidly it effectively sees the radially symmetric 
finite radius Fourier modes as the minimizing modes. 

C. Zl and Z2 potentials 

Lately, two new potentials for monatomic glass like 
formers in the bulk were found [20]. In their presence, 
the atoms display non-compact arrangements of 13-atom 
icosahedra and a hint of fragile glassy dynamics in simu- 
lations of the cooled liquid state. These are [20] 

V{\r\) = aexp[-ar]/r^ cos{2kFr) + h{<j/rY + Vo, (31) 

with certain numerical values for the various parame- 
ters {a, a, /cf, cr, n, Vo} for r > Vc- This model can be 
related to a superposition of Kac-type potentials. As 
in all glass models that we investigate, the minimizing 
Fourier modes of v{k) occur on shells. Like other po- 
tentials, V(r) naturally leads to icosahedral like ordering 
over large non-compact regions. In the continuum limit, 
these potentials are akin to non-Abelian backgrounds. 
The extended space that they occupy offers a genuine 
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FIG. 5. The minimum energy clusters of tlie Zl minima. 
Tiie number next to eacli cluster denotes tlie number of atoms 
in it. Tliis figure is reproduced from Doye et al. 

example (not merely that of small 13 particle clusters) 
of relatively large systems (~ 10^) with geometric frus- 
tration. The local maxima of the potential inhibit closed 
packed structures and lead to polytetrahedral structures 
which may be divided into tetrahedra with atoms at their 
corners. Local icosahedral order involves no disclinations 
but as the clusters enlarge, negative disclinations must be 
introduced. The disclination lines cannot end abruptly 
and must form closed loops or traverse the cluster. As in 
Frank-Kasper phases, the packing can be described by a 
network of disclination lines threading an icosahedrally 
coordinated medium. Clusters of various sizes are de- 
picted in Fig. (5) from [20]. For the two potentials for the 
largest clusters, networks and chains of icosahedra are re- 
spectively preferred. Upon cooling, in silico, a liquid of 
^ 10'' particles, extended icosahedral clusters form with 
long time for local rearrangements. These potentials pro- 
vide a natural realization of the notions underlying geo- 
metric frustration (here with icosahedral like order over 
a large range), and a naturally related glassy dynamics 
(slow rearrangement of icosahedral aggregates inhibiting 
crystallization leading to better glass- formers). 

IX. CORRELATORS AND DISORDER LINES 

A. The n ^ oo limit 

The correlation functions (and much of the physics) 
are governed by the "thermal" dynamics of the zeros 
and branch cuts of the energy spectrum in the complex 
k plane. Generically, in a system with many compet- 
ing interactions- with either direct competition (as in the 
Coulomb frustrated ferromagnet) or indirect (e.g. non- 
Abelian backgrounds)- the Green functions will display 
numerous poles and/or branch points. These singular- 
ities, generally, adiabatically shift their location in the 
complex plane as a function of temperature (the chemi- 
cal potential fi or mass gap r). Exactly how they evolve 
in the complex k plane determines, to a large extent, the 



nature of the transitions and crossovers that occur. 

In all of the models that we consider- all models with 
a minimum of the kernel on a shell in |fc| space- we find 
that at low temperatures the branch point or pole of rel- 
evance has a finite real part leading to oscillations. This 
is expected as at low temperatures, the relevant Fourier 
modes are on the shell |fc| = g and the system must dis- 
play oscillations. As temperature is increased, the mod- 
ulation length becomes longer and longer. The evolu- 
tion of the modulation length can be either continuous 
or discontinuous. For systems with finite range interac- 
tions (including non-Abelian backgrounds), the number 
of poles and branch-points is conserved as a function of 
temperature- no two (or more) singularities merge in the 
complex plane. In such systems, the modulation length 
may, at times, jump from one value to another. This 
occurs when one pole or branch point approaches the 
real axis more than another singularity (and consequently 
has a smaller imaginary part). The relevant asymptotic 
modulation length then jumps to the singularity which 
is closer to the real axis. In our toy S0(3) actions, 
this indeed happens: At high temperatures, the sys- 
tem is, asymptotically unmodulated. Below a crossover 
temperature, the system jumps and starts to exhibit fi- 
nite modulations associated with singularities of another 
eigenvalue. Naively, one might expect the branch points 
present in the eigenvalue spectrum to immediately lead 
to algebraic (i.e. non-exponential) decay of correlations 
at all temperatures. This, however, is not what occurs. 
Cancelations arise and correlations decay exponentially 
at all finite temperatures (as expected). 

In systems with competing interactions of different 
longer range interactions, the poles merge frequently. For 
instance, in the Coulomb Frustrated Ferromagnet, the 
poles move along the imaginary axis at high tempera- 
tures. At a certain crossover temperature, they merge 
and hit the "unit" circle. From there, the poles bifur- 
cate in pairs to form two conjugate pairs of poles on the 
unit circle leading to finite modulations which become 
of shorter and shorter length (and ever larger correlation 
length) as the temperature is decreased. The merger of 
the poles at the crossover temperature leads to a "disor- 
der line" like effect. The thermodynamic functions are 
analytic at this temperature yet their explicit real func- 
tional form changes (just an analytic continuation from 
one natural set of functions to another) . In systems with 
long range interactions (e.g. the Coulomb Frustrated Fer- 
romagnet of Eqn.(27))), the modulation lengths tend to 
decrease as temperature is decreased; the converse oc- 
curs for the canonical interactions (fc^) thwarted or frus- 
trated by an even shorter range interactions (e.g. k'^). 
In physical terms, the crossover (disorder line) tempera- 
ture is related to a temperature that marks the gradual 
crossover from a one dimensional system (the relevant 
density of states at low temperatures is indeed one di- 
mensional) to a canonical three dimensional system at 
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high temperatures. For exphcit expressions for the cor- 
relation and modulation length crossovers the reader us 
referred to Appendix(C). A discussion of the thermody- 
namic crossovers appears in Appendix(D). 

B. A Josephson like length 

At very low temperatures the real space correlation 
function G(x) decays differently in different regimes. A 
0{l/n) calculation reveals an algebraic cross-over similar 
to an inverted Josephson like regime. Specifically, the 
pair correlator 

G{x) ~ |f for^j «: \x\ < -Ld < ^, 

G{x) ~ l^r^ for \x\ «; Ij, (32) 

with the Josephson like length 

ij = n/(2E(0)). (33) 

Evaluating matters exphcitly ( [11], [24]), we find that 
the Josephson like correlation length £j scales as the in- 
verse of the correlation length in the frustrated system. 
Here S(o) is the 0(n°) self- energy, ^ is the correlation 
length, and Ljj, the domain size, is {2Tr/\Ke{k-}\), with 
"Re" denoting the real part. More generally, within the 
general SO (4) theory, k- is replaced by the root of the 
eigenvalue corresponding the largest correlation length 
(smallest imaginary part). This length scale appears 
for all actions in which the relevant interaction kernel 
is peaked about a (d-1) dimensional manifold \k\ = q> 
about which they are quadratic. 

X. A GENERALIZED MERMIN- WAGNER 
INEQUALITY 

Here we report on an extension of the Mermin- Wagner 
theorem which naturally links the magnitude of the or- 
der parameter to constraints on the characteristic deco- 
herence or relaxation time. We generalized the Mermin- 
Wagner inequality to all two dimensional systems with 
an analytic interaction kernel and to a host of three di- 
mensional systems in which a certain integral diverges. 
The inequality that we derive reads 



2|m,|2T(rg) < 1 



with rrig the magnetization at the minimizing momen- 
tum 5, and (with ?i = 1) r^: is an integral having the 
dimensions of time reflecting the characteristic decoher- 
ence or relaxation time in a quantum system suffering 
a boost of momentum k (scaled magnon liftime). If the 
characteristic relaxation time diverges the system does 



not admit, at any finite temperature, a fluctuation spec- 
trum about an ordered state. As we will later show be- 
comes glassy (its configurational entropy diverges). The 
characteristic decoherence time 

l/{Af E)^rj:, (34) 

with ' E the change of energy of the system due to a 

boost of momentum k or scaled spin- wave energy. When- 
ever the integral 



(27r)<«y v{k + q) + v{k-q)-2v{q) 

diverges then we may not obtain the low temperature 
T = 0"*" magnon dispersion by assuming a nearly per- 
fectly order state. In the above, g G M is any vector 
on the minimizing manifold (M) in momentum space. 
For further details the reader is referred to the Appen- 
dices. In Appendix(E) we generalize the Mermin- Wagner 
inequality to all rotationally invariant two dimensional 
system with an analytic interaction kernel in momen- 
tum space. In Appendix(F), we derive inequalities in 
three dimensions and interpret the extended Mermin- 
Wagner theorem in terms of relaxation times. When 
multiplied by the temperature, the integral appearing in 
the generalized Mermin- Wagner inequality (Eqn.(35)) is 
a strict lower bound on thermal fluctuations. Spin sys- 
tems with a non-trivial minimizing manifold in k space, 
much like spin ladders, display an interesting even-odd 
effect. When we analyze the fluctiiations of spins having 
an even number of components (n) , we find precisely the 
integral of Eqn.(35) appearing as the relevant thermal 
fluctuation integral. For a system of spins having an odd 
number of components, we find a more divergent (for on- 
shell minima in Fourier space) fluctuation integral. For 
further details the reader is referred to [33]. The bound 
that we derive by a generalized Mermin- Wagner inequal- 
ity is a strict lower bound on the fluctuations. When 
the number of spin components (n) is odd, the thermal 
fluctuations are much larger, and ordering is even more 
inhibited. 

XI. ENTROPY AND DEGENERACIES 

The generalized Mermin- Wagner inequality suggests, 
divergent decoherence times are tied to a vanishing order 
parameter. Physically both effects stem from one origin- 
the complicated energy landscape that these systems pos- 
sess. The high degeneracy and near degeneracy of low 
lying states makes these systems vulnerable to thermal 
fluctuations. Metastability is at the heart of the matter. 
It is not surprising that the configurational entropy may 
be extensive at high enough temperatures. 
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A. Ground State Degeneracies 

The ground state degeneracies of this system are ob- 
vious. As the enc'rgy dispersion in these systems is ef- 
fectively one dimensional (it depends only on the radial 
coordinate \k\), whereas the system is d dimensional, the 
ground state degeneracy should be exponential in L'^~^, 
with L the system size. These statements are exact in the 
spherical (large n) limit where every normalized superpo- 
sition of minimizing modes constitutes a ground state. As 
the number of minimizing modes is 0(L'^~^) [33], the pre- 
vious exponential degeneracy follows. For special Ising 
models, we can rigorously prove [33] that the degener- 
acy is bounded from below by 2l^l where \M\ denotes 
the number of points on the minimizing manifold which 
is characteristically 0{L'^~^). For a similar analysis for 
0{n) groimd states with a trivial bound see [33]. 

In the largo n limit all of this becomes exact and the en- 
tropy scales as the surface area leading to a holographic 
like ground state degeneracy. 



B. Replica Calculations 

For metastablc states at higher temperatures the prob- 
lem is harder and we were not able to derive rigorous 
results. A replica calculation within the self-consistent 
screening approximation reveals that the configurational 
entropy in all of these systems: continuum non-abelian 
gauge backgrounds, Kac-type interactions [22] , Coulomb 
frustrations [21], and in general any system in which 
the interactions are peaked on (d-1) dimensional shell 
in Fourier space, is extensive. 

The case of Coulomb frustration was investigated by 
[21] and generalized to Kac-type step- function interac- 
tions in [22]. Here we note that the minimal ingredient 
in all of these calculations, the requirement that the un- 
perturbed Green's function be of a simple Lorentzian or 
delta- function form at low- temperatures [22,24] 



Go(fc) ^ 



+ (A;2 - ^2)2 ■ 



(36) 



peaked about the manifold of minimizing modes = g 
(Fig.(l)). As in [22], here Z denotes the weight of the 
peak and ^ the correlation length or inverse mass gap. 
For any unperturbed Green's function which may be 
written as a Laurent series in fc^, we may approximate 
the Green's function in the vicinity of its maximum by 
Eqn.(36). Eqn.(36) is the form analyzed long ago by 
Alexander and McTague [34] in the context of the liq- 
uid to solid transition. The bare Green's function Gq 
is the approximate structure factor for a liquid with / 
the average inter-atomic distance. If no cubic terms fa- 
voring crystallization are present in the action (i.e. if 
the transition to the crystal is "avoided"), and perform a 



self-consistent screening approximation (as done by [21]), 
then we will find a glass at low temperatures. An evalu- 
ation of the configurational entropy [22], [21] 



-I 



(27r)3 



S[ 



(uT)-i + G^G 



]} (37) 



with S[x] = —X — ln(l — x) shows that the entropy is 
extensive within a broad temperature range (T^ < T < 
Ta = Tc{q = 0)/(q'7r^ -|- 1) [21] and a vanishing configura- 
tional entropy density Sc{Tx)/V = defines the Kauz- 
mann temperature T^). In the above, F is the Edwards- 
Anderson order parameter limt^ao {4'{~ t)4>{k, 0)) , G is 
the standard correlation function, and denotes a con- 
volution in Fourier space [22]. 

We now envision inverting the problem. Instead of 
starting from a certain fixed model and working our way 
through to find the configurational entropy, we go back 
to ask what sort of systems will exhibit extensive con- 
figurational entropy (i.e. an exponentially large num- 
ber of metastable states) by the via provided by calcu- 
lations that invoke the approximation of Eqn.(36). The 
answer is that the "universality class of systems" with 
Fourier transformed kernels v{k) having their minima on 
the shells |fc| = g > 0) and which are analytic in the 
vicinity of these minima, can have their Green's functions 
approximated by Eqn. (36) (e.g. non-Abelian background 
fields, myriad translationally invariant frustrated models 
with rotational symmetry). Once we assume this form, 
we may replicate the replica calculations of [22], [21] word 
for word to find that the configurational entropy is pro- 
portional to the volume Sc a V, leading us to conclude 
that the system displays an exponentially large number 
of metastable states. The calculations of [21], [22] can 
be reproduced with slightly modified expressions for the 
onset temperature Ta of glassy dynamics in various sys- 
tems. Within all of these systems, including the domi- 
nant contribution from the lowest eigenvalue appearing 
in non-abelian gauge theories of glasses earlier, we find 
identical results. Repeating this calculations for temper- 
atures T <Ta, the self-consistency conditions detailed in 
[21] are found to allow a nontrivial glassy state wherein 
defects may not, in all cases, venture more than 0{q~^) 
at arbitrarily long times. 

For the SO (4) action of Eqn. (7), the lowest rank n = 12 
representation has N = 338 individual real components. 
Within the large N ^ 1 limit, each individual complex 
eigenvector has a correlation length ^ associated with 
it and the configurational entropy 



_ 'y ^ ^M,i,n 

MA,n 



(38) 



where 5^'*'", the configurational entropy associated with 
each field ^ {x) (the Fourier transforms of the eigen- 

mode amplitude ct^ i{k)) is given by the self consistent 
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screening approximation [21], [22] for a single scalar field 
(and carried in works till now by a large N calculation, 
setting TV = 1 at its end). The field associated with the 
lowest eigenvalue, (pmin: has its minimizing wave- numbers 
on a finite radius, q > 0, shell and its configurational 
5™^" is extensive within the regime T^^" > T > T™^". 
Trivially generalizing earlier calculations [21], the config- 
urational entropy originating from the lowest eigenvalue. 



o3. 

4-77 



V 



-(1 



^min \2 



2 / ._. ^min 

— u z — 



)2+ln(l-(l 



^min 



(39) 



Here [21], [22] Ki„i„ satisfies 



8e2 



(1 



i/ /^min) 



i-(i. 



n/'^min)^ 



with e^in = ^mfn/«'4in and V the volume. ^^'"(Tc) = 
C^«min' with the constant C ~ 1.18 x lO'^. T^'^ may 
be for the single scalar field (^min [21], [22]. As the con- 
figurational entropy is a sum over all eigenvalues, Sc is 
generally extensive over a larger range than T™'^ suggest. 



XII. SPATIALLY NON-UNIFORM SYSTEMS 

Thus far, wc discussed ideal uniform frustration. We 
now illustrate that some of our ideas and calculations 
might not be so restricted. Any kernel Vij{x,y) in the 
general two body Hamiltonian of Eqn.(3), becomes diag- 
onal by some unitary transformation. The Fourier modes 
are the eigenmodes of V in the ideal translationally in- 
variant setting. In general, Vij {x, y) becomes diagonal in 
another complete orthogonal basis ({|u)}): 

{Ua\Vij\U0) = 5al3{Ua\Vij\Uf)) = SapVij{Ua). (40) 

We may now examine the minimizing manifold in u space 
for the lowest energy eigenvector. Repeating calcula- 
tions, we find that if this surface if {d — 1) dimensional 
(e.g. a shell) and the relevant Vij{u) is quadratic in its 
environs then generally, for large n, Tc = 0. The rig- 
orous, large n, holographic ground state degeneracies of 
the previous section and the large n configurational en- 
tropy analysis retain the same character. In many of 
the expressions above, the momentum index k is merely 
replaced by the more general diagonal basis index u. 



scale. As a result, it hard by principle to make precise 
empirical statements in the vicinity of the "glass transi- 
tion" itself (where the relaxation times, according to com- 
mon wisdom, becomes infinite). We now discuss possible 
non-rigorous derivations for both the Vogel-Fulcher (VF) 
form (subsection A) and the avoided critical fit (subsec- 
tion B). These two fits are different. There currently is 
no sufficient experimental or mathematical proof of any 
these forms. It is, ab initio, possible that both forms and 
their some of their assumed underlying theories are two 
aspects of the same problem. For instance, it is possi- 
ble that the entropic-droplet/VF description is valid at 
very low temperatures whereas the avoided critical ex- 
pression has its justification at higher temperatures (not 
far from the avoided critical temperature itself). Our cal- 
culations indicate that the theories of geometric frustra- 
tion in glasses display both an avoided critical point and 
extensive configurational entropy. Albeit their suggestive 
character towards avoided critical and/or VF dynamics, 
none of these fits rigorously follows from the thermody- 
namics that we investigated. 



A. A Possible Derivation of Vogel-Fulcher Dynamics 
In Non-Abelian gauge backgrounds and others 

The basic premise of the VF fit is that the relaxation 
times seem to diverge at a temperature which correlates 
well with the intercept (Tr-) of the extrapolated entropy 
of the supercooled liquid with that of the solid. As out- 
lined in the previous section, as the system is cooled 
down, the configurational entropy first becomes exten- 
sive Ta (at which free energy minima first appear); at 
Tk these minima become stable and the configurational 
entropy vanishes. Once we find an exponentially large 
number of metastable states for Ta > T > Tk in non- 
Abelian backgrounds and other systems (by the replica 
calculations of the previous section), we may invoke the 
analysis of c;ntropic droplets [25], [21]. This leads to char- 
acteristic energy barriers 



a {TSc 



(41) 



Linearizing the extensive configurational entropy, 

Sa{T ^ T+) ~ V{T/Tk - 1), (42) 

this now leads to VF dynamics [25]. In the VF fit, the 
relaxation times 



XIII. GLASSY DYNAMICS 

In glasses, the relaxation times may increase by 15 or- 
ders of magnitude over a temperature range of lOOK. It is 
impossible to obtain data below the glass transition tem- 
perature. All data has to be analyzed on a logarithmic 



t^cxp[DTk/{T-Tk)]. (43) 

Here, Tk is the Kauzmann temperature at (or above) 
which an "ideal glass transition" would occur. At Tk 
the extrapolated entropy of the liquid undergoes a crisis. 

Thus, once we show that the lowest eigenvalue of the 
non-abelian action has a minimum on a shell then we 
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may non-rigorously derive VF dynamics. According to 
this line of logic, many of the translationally invariant 
systems with minimizing momenta \k\ = q > 0: 

• Rigorously, do not order (at least in the usual sense) 
at finite temperature (the transition is narrowly avoided 
for 1 > > 0). 

• May display many of the characteristic thermody- 
namic signatures of narrowly avoided phase transition 
(e.g. a crossover of the thermodynamic functions at tem- 
peratures "Ti" which merge with the critical temperature 
of the unfrustrated system in the limit of zero frustration 
(i.e. when q ^ 0)). 

• Rigorously, have numerous degenerate ground states, 
and, on a less rigorous footing, by an insertion of Eqn.(36) 
into Eqn.(37), have macroscopic configurational entropy 
(an exponentially large number of mctastable states). 

• The extensive configurational entropy of geomet- 
rically frustrated systems may lead to glassy (Vogel- 
Fulcher like) dynamics at low temperatures. We may now 
non-rigorously derive VF dynamics by fusing standard 
approximations [25] with the extensive configurational 
entropy results of the replica calculations (Section(XI)). 

B. A Derivation of an Avoided Critical Point in the 
Avoided Critical Dynamics Ansatz 



component theories of glasses {n > 12 in the (n + 1)^ 
complex components of the icosahedral multipole mo- 
ments Qn)- In Fig- (2), the ideal freezing temperature 
Tif = Tc{l/R = 0) is the transition temperature on the 
ideal template that is narrowly avoided in flat Euclidean 
space. This provides serious conformation to the notions 
imderlying the avoided critical fit. We further proved, 
by a generalized Mermin- Wagner theorem, that no stan- 
dard dispersion relations about nearly ordered states can 
exist at arbitrarily low temperatures for a host of sys- 
tems, all of which display a minimum of the interaction 
kernel at finite wave-numbers compromising a high di- 
mensional manifold in k space. Lately, numerical results 
[27] for clock realizations of a lattice version Hamiltonian 
of Eqn.(27), indeed show agreement with the form pre- 
dicted for the scaling form of the free energy barrier [12] 
for this Coulomb frustrated model (Eqn.(44)). Further- 
more, as we will now show, if we compare the magnitude 
of the avoided critical temperature {Ti f ) as inferred from 
the relaxational dynamics fit the free energy barrier, we 
will find that it agrees very well with what we would ex- 
pect on simple energetic grounds. We provide further 
empirical impetus for assumption (2). 

C. Other Relaxational Forms 



As heralded in [12], perhaps nature follows another ap- 
proach. One of the troubling (or surprising) aspects of 
the "glass transition" (if it is indeed a thermodynamic 
transition) is that although the dynamics undergoes a 
stupendous change- it may increase by orders of magni- 
tude over a small temperature window- the smooth ther- 
modynamic functions do not seem to know about any 
impinging transition. My coworkers and I [12] raised the 
possibility that the dynamics is governed by an avoided 
critical point. This viewpoint has two main assumptions: 

1) The system would have undergone a phase transi- 
tion in high enough dimensions but due to the frustration 
encountered in low dimensional Euclidean space it can- 
not. The phase transition is narrowly avoided leading to 
the phase diagram of Fig. (2) . 

2) The dynamics is tied at its core to the thermody- 
namics. The exponential increase in the relaxation times 
as the "glass transition" is approached is simply a reflec- 
tion of a super-Arrhenius type dynamics r ^ exp(£'/T) 
with the free energy barrier E a function of the re- 
duced temperature as measured relative to the; avoided 
critical crystallization temperature on the unfrustrated 
high dimensional template (e.g. the 4-sphere on which 
icosahedral ordering can occur in the case of the simple 
monatomic system). 

Here, we proved that in the large n limit, assump- 
tion (1) for the geometrically frustrated actions intro- 
duced for glasses by many [1]. Our calculations yield 
the phase diagram of Fig. (2) for the N > 338 » 1 real 



Following this work, based on a new Monte Carlo anal- 
ysis, [28] suggested that scalar systems with an on-shell 
momenta minimum display a pronounced slowing down 
yet not nearly as dramatic as that in true glasses. The 
dynamics is consistent with Hartree results. If true, when 
fused with the contents of this paper, this finding suggests 
that non-Abelian gauge background trigger a novel slug- 
gish dynamics on their own right. This new dynamics, 
albeit very interesting, will have little to do with actual 
structural glasses. As detailed in later sections, strTic- 
tural glass theories must account for dynamics borne by 
topological defects (ignored in the uniform gauge back- 
ground treatment). Within the relevant range of inves- 
tigated parameters [29] , a critical analysis of the data of 
[28] does not rule out glassy dynamics. It is possible that 
none of the forms presented (including the Hartree form) 
is correct for the continuum models of earlier sections. 
Nevertheless, all current simulations [27], [28] unambigu- 
ously point at sluggish dynamics implying very slow dy- 
namics in non-Abelian backgrounds. Ingoring topolgo- 
ical defects all Landau- Ginzburg theories of liquids and 
crystals have a quadratic term with v{k) having minima 
on the shell = g (augmented by cubic terms in crys- 
tals and only quartic and higher in liquids) [34]. This 
suggests that if explicit topological d(^f(x;ts are ignored 
then the most natural continuum actions for supercooled 
liquids are those which we investigate. 
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(1-T/t") 



FIG. 6. This figure is from D. Kivelson et al. Here we 

fitted the relaxation times of all known glass formers to 
the super- Arrhenius form of Eqn.(44). The horizontal axis 
(1 — T /T^^) gauges the reduced temperature as measured rel- 
ative to the transition temperature T^^ on the ideal curved 
template. 

XIV. PREDICTING Tif IN THE AVOIDED 
CRITICAL FIT FOR GLASSES 

We now examine the avoided critical dynamical fit 
of glasses [12] and present a comparison between the 
average avoided critical temperatures in these fits and 
the actual expected from simple energetic considerations 
(which seem to coincide remarkably well). 

We saw that the simplest gauge theories of metal- 
lic glasses display an abrupt drop in their transition 
temperatures. The temperature of crystallization in 
3 — dimensional Euclidean space is depressed well below 
the crystallization temperature on the curved ideal tem- 
plate. This discontinuity might be of profound impor- 
tance: the thermodynamics (and consequently the dy- 
namics) might be governed by the proximity to the "ideal 
freezing temperature" T*-^- the temperature of crystalliza- 
tion on the ideal curved space template. In [12], we 
demonstrated that the following empirical fit works re- 
markably well for all known glass formers: 

r ~ Too exp[£(T)/(/esT)]; E{T) = + {atf/^0{t) (44) 

Here, t = [T^f - T)/T^f , the relaxation time r is de- 
duced from viscosity and time dependent heat capacity 
measurements, and 9{t) is the Heavyside theta function. 

We may expect ksT^^ to be indicative of the cohesive 
energy of the ideal curved space crystal, and the natu- 
ral melting temperature scale (fesTnie;*) to be correlated 
with the crystallization energy of the real (space) struc- 
ture. Obviously entropic efi^ects are important (e.g. those 
present in the BCC — > FCC transition [34]). Amongst 



these liquids, the average value of a; = [T*-^ — Tm]/T„i 
is 9.59%. We may compare this to the value com- 
puted by Frank [26] for the relative cohesive energy, 

l^^tcosahderai_^FCCy^FCC^ foj. ^ simple Lcnuard Jones 

liquid, which leads us to expect (x) = 8.4%. The prox- 
imity of these two values is misleading, x varies dramat- 
ically from one liquid to another- the ideal crystalline 
structures, and their respective energies change radically 
from one material to another. Cautious optimism is still 
called for. If the theory were reliable and the parame- 
ter T*-^ obtained from the fit could indeed be interpreted 
as the freezing temperature on the ideal template, one 
would expect that, for all liquids, a; > 0, as the structure 
on the ideal template (where the ideal minimum energy 
local structure could be extended in all directions) is by 
its very definition more stable than (or at least as stable 
as) the corresponding structure constrained to reside in 
Euclidean space. 



glass former 


X = [T^f - Tm\/Tm 


n-butyl benzene 


7.02% 


triphenyl phosphite 


8.2 % 


isopropyl benzene 


20.6 % 


propylene carbonate 


7.73 % 


salol 


-4.4 % 


dibutyl phthalatc 


21 % 


o-terphenyl 


5.7 % 


s-trinaphthyl benzene 


9.1 % 


n-propanol 


30.6% 


a-phenyl-cresol 


-10.3 % 


glycerol 


9.9 % 


glass former average 


9.59 % 


simple theory 


8.4% 



XV. AVOIDED CRITICAL DYNAMICS AND 
MICROSCOPICS 

Ideal order cannot tessellate Euclidean space and it 
may do so only on a higher dimensional surface with en- 
ergetics that seem to concur with the picture that we 
have in mind. In order to make progress toward un- 
derstanding the dynamics we need to explicitly sec how 
things transpire in real Euclidean space. We may adopt 
two trivial points of view: either the action is pinned to 
be that of a system in a fixed non-Ahelian backgroimd, 
or the distribution of matter alters to the geometry by 
virtue of the topological defects that it carries in which 
case we have a full blown matter coupled gauge theorjj- an 
elastic theory with disclinations where the disclinations 
are not pinned down. The ideal order on the curved high 
dimensional template needs to be mapped into fiat space 
by introducing the disclinations of Rivier to compensate 
for the lack of curvature in Euclidean space. 
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A. Fluctuating Geometry- 
Let us assume that the ideal curved order is neutral- 
ized in flat Euclidean space by the insertion of a uniform 
disclination density (giving rise to an effective negative 
curvature to compensate for the positive curvature on the 
torn pieces of the ideal sphere). When glued together the 
curved pieces of the fractured sphere form disclinations 
at their common junctions. This is much like the seams 
that we encounter in orange peals that need to glued 
together somehow to form a flat surface. The disclina- 
tion lines can easily be traced by looking at the locus 
of all points that have a lower or higher than usual co- 
ordination number as compared to that present in the 
ideal icosahedral tessellation. Formally, they can be seen 
by decurving processes relating the curved template to 
flat space. To a zeroth order approximation, the addi- 
tional energy of the decurved structure (an elastic strain 
energy) amounts to the total length of the disclination 
network. In the context of decurving the {3,3,5} poly- 
tope, the net number of disclinations is dictated by the 
number of sites having a coordination number different 
from the ideal (12-fold coordinated vertices). An anal- 
ysis along this counting may be found in [1]. If we let 
the field (j) denote the density of disclination lines then 
in our analysis such a zeroth order effect amounts to a 
chemical potential term with the chemical potential 
/X now denoting the energy per unit disclination line. 

Disclinations (or disclination pairs- dislocations) inter- 
act by direction dependent Coulombic like forces. In a 
two dimensional XY model like setting (e.g. vortices in a 
thin film superconductor) where the disclinations amount 
to additional twists that the phase must undergo when 
parallel transported about their cores. Such fractional 
vortices interact Coulombically. Two fractional vortices 
of opposite charges undo the damage that the other does 
at long distances. In three dimensions, the energy would 
be minimal for such wedge disclinations (or fractional 
vortices) lying in the same plane. In the tensorial elastic 
setting, such angle dependent forces undergo only simple 
detailed change. The Blin interaction between disloca- 
tions is a decorated Biot Savart interaction [32]. 

The same occiirs in gravity. In 2+1 dimensional grav- 
ity, masses amount to conical defects which cut and re- 
move (or insert) additional angles in the plane much like 
disclinations. If we disregard "unimportant" details, try- 
ing to get to the heart of the matter by noting that the de- 
fect densities interact Coulombically, then we may imme- 
diately provide a possible derivation of the Coulomb frus- 
trated Hamiltonian of Eqn.(27). We now try to give the 
reader a quick intuitive grasp of the microscopic deriva- 
tion that what lacking till now (as still is). Let (f) denote 
the defect density. In reality </) is a tensor by let's disre- 
gard details and merely see how things may work. For 
simplicity, we set the field ^ be the local scalar curvature. 



It is positive inside the broken fragments of the ideal 
sphere that have been fiattened out; the net curvature 
is superposed on a neutralizing background that ensures 
that the net curvature is zero (as in the Frank-Kaspcr 
phases). Once we make the correspondence among the 
ideal sphere fragment densities and the disclination densi- 
ties, eqn.(27) may, trivially, follow. The disclinations (or 
disclination pairs) interact Coulombically and the density 
toward ideal (high dimensional spherical) order within 
each resulting spherical fragment is captured by the gra- 
dient term. These two terms compete leading to exten- 
sive entropy at low temperatures and glassy dynamics. 
However, the dynamics is not trivially of a Vogel-Fulcher 
form. 

To give the discussion a more detailed flavor we now 

give some of the standard expressions. Denoting by u 
the displacements of atoms from their ideal equilibrium 
positions, the dislocation density reads 

aij = CikidkdiUj. (45) 

In a gravitational language, the dislocation density is, in 
simple disguise, the connection determining the change 
the of a vector due to parallel transport alone, 

Oiij = CikiTkij (46) 

with the connections 

Lijfe = didjUk- (47) 

The connection between disclinations and gravitational 
curvature is even more transparent. With 

^ij = ^{9iUj - djUi), (48) 
the disclination density is 

&ij = ^eikiejqsdkdiuiqs, (49) 
while the curvature tensor is 

Rijki = {didj - djdi)wki. (50) 

An explicit calculation for a wedge disclination corre- 
sponding to the insertion of an additional angle (of mag- 
nitude a) in between the lips of a cut in the 12 plane 
shows that 

i?i2i2 = -aS^ix^) (51) 

where d'^{x^) denotes the delta function in the plane 
for the disclination line x''^ threading three dimensional 
space. Thus, the curvature tracks the disclination lines 
and its magnitude is proportional to the disclination an- 
gle. The disclination lines of Rivier [15] become mass 
sources. For a positive inserted angle, the scalar curva- 
ture is negative. If we wish to squash the angles to close 
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the gaps such as those appearing when trying to piece 
together five tetrahedrons around a common edge, we 
insert a negative angle. The curvature, the full contrac- 
tion of the curvature tensor, 

R = R';Z (52) 

is positive on a sphere (scaling as the inverse squared 
radius). The Blin interaction between dislocations 

E=^Jd^x (fx' aij{x)[K]avy, (53) 

with the distance dependent part of the kernel [K\ 

-{fE]^l22ll) + _A_eyj.e;,^.,^,aj.a;.,r, (54) 

where ji and v are the elastic constants, and a the dis- 
location densities of Eqn.(45). All interactions between 
elastic defects arc of a Coulombic character which has 
been tensorially flavored up. We can easily see why this 
must be so: Elasticity is centered around a second order 
gradient expansion, 

1 

^ii — + djUi), (55) 

with u the displacements about an equilibrium position. 
When Fourier transformed and inverted, the Green's 
functions through which the topological defects commu- 
nicate always give rise to ~ j.'^-d interactions. In 
two dimensions, the interaction between the dislocations 
is, not too surprisingly, logarithmic. As all derivations 
that we have been able to devise (apart from the gener- 
alized Mermin- Wagner theorem and the construction of 
0(n) ground states) have been approximate, we examine 
numerical results. Lately, computations [27] for clock re- 
alizations of a lattice version Hamiltonian of Eqn.(27), in- 
deed show agreement with the form predicted for the scal- 
ing form of the free energy barrier [12] for this Coulomb 
model (Eqn.(44)). We have given simple arguments for 
the form of this fit, although a real theory is still lack- 
ing. As noted, the geometrical frustration breaks the 
glass into clusters. In an earlier publication [12], we 
have found that assuming the usual Ornstein-Zernicke 
form for T < T^, within our "standard model" , cumulant 
crossing occurs when the size of the domains is larger 
than Rjj — '^ith the correlation length of the 

unfrustrated system which diverges at the (avoided) crit- 
ical temperature T^. A somewhat similar result (with 
^0 replaced by the correlation length ^ of the frustrated 
system) followed from our 1/n calculation [11] for the 
Josephson like correlation length. This 1/n calculation 
result undergoes no change for any of the models that 
have their free energy minima on a shell in Fourier space 
which becomes sharp at low temperatures. 



To shortly examine viable extensions later on, we 
sketch our derivation in [12]. Usually, the specific heat 
exponent a is small and consequently from the hyper- 
scaling relation: 

vd = 2-a, (56) 

we expect i/ ~ 2/d. The free energy barriers that govern 
the low temperature dynamics involve the flipping of do- 
mains (of size Rd)- We assume that all is governed by 
the relative proximity to the avoided critical point. As 
in Eqn.(44), we consider temperatures below the avoided 
critical temperature T < T*-^ and denote by t the ab- 
solute value of the relative reduced temperature. The 
correlation length in the unfrustrated system S,o ~ t~'^. 
The surface tension in the unfrustrated system 

so 

with B a constant. The free energy barriers that we 
might find by breaking the medium into unfrustrated 
patches of size Rjj (the geometrical frustration in this 
approximation merely breaks the unfrustrated medium 
into finite patches (domains) of scale Rd), 

Ec^aR^-'c^BTifi^f-^ (58) 

Employing our estimates [11,12] Rd ~ ^o"^ and i/ ~ 2/d, 
E ~ t4(d-i)/d_ ^59^ 

Substituting = 3 this gives an exponent of 8/3 which 
which works best in our fit. It should be noted, however, 
that if we assume that 

E{T) ^E^ + {atf (60) 

then for we can obtain reasonable fits with 7/3 < a; < 3. 

B. Fixed Non-Abelian Background 

Here no derivation of Eqn.(44) seems possible along 
the lines of the previous derivation. We may read off the 
dispersion from the lowest eigenvalue. For instance, in 
our toy model for the S0(3) action, A_(fc) = //, + k'^ + 
36*2/2 - \\/TQkWT^ is the relevant dispersion. For 
k ^ 9, (distances small compared to the zero tempera- 
ture modulation length- the radius of the ideal sphere) 
the square root may be expanded and to lowest order, the 
dispersion is trivially {k"^ — 2k6 -\- n) . In a spin model, such 
a dispersion corresponds to a short range nearest neigh- 
bor ferromagnetic interaction [k^) frustrated by a longer 
range anti- ferromagnetic interaction {—29\k\). Unfortu- 
nately, the frustrating interaction is not of a much longer 
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algebraic decay than the unperturbed ferromagnetic in- 
teraction. In a cumulant expansion similar to the one 
that we performed in [12], we now find that the charac- 
teristic domain size is Rd ^ When inserted back into 
our derivation for the free energy barriers this gives phys- 
ically impossible dynamics: E ~ i2(i-d)/d^ £j.gg energy 
barrier that quickly diminishes as we go away from the 
avoided critical point. For k <^ 9 this worsens. We can- 
not see a natural generalization of the frustrated domains 
of the previous section which provides sensible results. 

To summarize, for non-Abelian gauge backgrounds, we 
can very easily derive Vogel-Fulcher dynamics yet we are 
unable to derive avoided critical dynamics of the form 
of Eqn.(44). Notwithstanding, these models do display 
avoided critical behavior insofar as their thermodynam- 
ics is concerned. Furthermore, at low temperatures only 
the Fourier modes close to the minimizing shell should be 
of any importance. If a low temperature fit of the form 
of Eqn.(44) works well for the Coulomb frustrated Ferro- 
magnet at small yet finite frustrations e, then Eqn.(44 
might also capture the low temperature dynamics of 
all models (including the fixed non-Abelian background) 
that have their minima on a shell of finite radius in k 
space. The approximate identification can be made by 
Taylor expanding the interaction kernels in the various 
cases to quadratic order and seeing that they all match. 

XVI. CONCLUSIONS 



5) We discussed how in a multitude of a systems (Non- 
Abelian backgrounds and all others with an "on-shell" 
minima in Fourier space), the derivable extensive config- 
urational entropy suggests Vogel-Fulchcr like dynamics. 

6) We illustrated how many of these notions can be 
extended to non-uniform system whereby the interactions 
are diagonal in a non-Fourier basis. 

7) We have shown how the dynamics may indeed corre- 
late with the avoided critical fit of [12] by comparing the 
melting temperatures with the fitted temperatures. The 
average values are in good agreement with our theory. 

8) We suggested a via for strengthening the link be- 
tween avoided critical dynamics [12] with microscopic 
consideration regarding interactions between disclina- 
tions. 
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APPENDIX A: N ^oo AVOIDED CRITICALITY 
IN THE QUANTUM ARENA 



To summarize, wc have shown that a multitude of 
translationally invariant systems having their minimiz- 
ing modes lie on a (d-1) dimensional shell in k space 

1) Do not allow for a magnon dispersion relation about 
an ordered state at any finite temperature whenever a 
certain (Mermin- Wagner like) integral is seen to diverge. 

2) We have noted that the generalized Mermin- Wagner 
inequality relates the magnetization (or its absence when- 
ever entropy dominates) and the average of the expected 
relaxation times within the system. 

3) The correlation and modulation lengths generally 
exhibit dynamics as a function of temperature (the imag- 
inary time width). Crossovers in the asymptotic form of 
the correlation functions do, in general, occur at a finite 
temperature which veers toward the critical temperature 
of the unfrustrated system in the limit of zero frustration 
(zero curvature in myriad geometric models). 

4) The systems exhibit an astonishing ground state 
and near ground state degeneracy. We may easily pro- 
vide rigorous lower bounds on the ground state degenera- 
cies and show, via replica calculations, that the configu- 
rational entropy is extensive. Physically, it is this mul- 
titude of degenerate and nearly degenerate states that 
makes the system so fragile to thermal perturbations. 
As a consequence, the ordering temperature the ordering 
temperature is very low. 



The single n — > oo parquet diagram (Fig. (4)) reads 

1 



nuT ^-v f (fik 



(Al) 



with the bosonic Matsubara frequencies ojm = '2nmT. 



where = y jj, + v{k). The sum can be evaluated by 
the method of residues which gives 

oo 

[to^-I-2/^]"^ = ^cothTry for y > 0. (A3) 



m— — OO 



Using this we obtain 



E° = i / -^^coth(% (A4) 



From the identity 



coth(-) = l + 2nB(Y) 



(A5) 



where nsiw^) represents the bosonic distribution function 
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ns(|) = [exp(|)-l]-^ 



Thus 



V.0 _ 



d'^k 2nB{ 



(fc)H 



(27r)" 



(A6) 



(A7) 



As /i ^ /imin = ^iKq), this integral diverges more 
strongly than in the classical case. This divergence re- 
mains down to zero temperature. The integrand is (up 
to a constant) the static correlation function. 



(fc)H 



v{k) + n 



(A8) 



APPENDIX B: THE SO(4) SOCCER-BALL 
MODEL ON A CUBIC LATTICE 

For 0(4) spins with 50(4) couplings in d = 3: 



(Bl) 



For four-component spins lying on a three-dimensional 
lattice (i.e. when the continuum limit is not taken) 
having 5*0(4) couplings "[— 5(;r) * S'(y)]" the matrix 
—Vij{k)/2 reads 



-E-' 

B -F 

O -G 

\e F G D 



where 

A = cos 6 cos ki + (cos k2 + cos ks) 
B = (cos ki + cos fca) + cos 6 cos k2 

2 

O = cos fc; + COS 9 cos fcs 

3 

D = cos ^ COS ki 

1=1 

E = isinO sin ki 
F = isinO sin ^2 
G = isin6'sinfc3 (B2) 

Note the appearance of a non-uniform sign (and non- 
trivial phases) linking the various spin components. 
These play the role of frustrating interactions present in 
scalar theories; in a high temperature expansion for the 
partition function Z, not all terms are the sum of dif- 
ferent contributions, all adding in unison with the same 
phase (and therefore also sign). The alternating signs 



of the resulting constituent terms are analogous to those 
generated in other frustrated systems by competing ferro- 
magnetic and antiferromagnetic interactions in the eval- 
uation of the partition function Z. For a nonzero value 
of the angle 6, Tc{6 7^ 0) = 0. The minimizing modes lie 
on a (d — 1) dimensional shell. 



APPENDIX C: EXPLICIT EXPRESSIONS FOR 
CORRELATION FUNCTION CROSSOVERS 

1. Uniform SO(3) and SO(4) backgrounds 

For the 50(3) action of [7] in d=2 dimensions, sum- 
ming over a set of (2^-|- 1) x (2^-|- 1) representations with 
corresponding rank spherical tensor order parameters: 



(CI) 



with {fiell } chemical potentials- Lagrange multipliers- set 
to secure overall normalization of {(f)i}. In fc-space, 

^ = ^ E E E \4{knKAm + M.]; 



with cl>em{k) = J2 "Kfc)eL(fc)- (C2) 



The self-consistent equation of constraint that the La- 
grange multipliers must satisfy is 



bT j 



(C3) 



The inverse critical temperature, , diverges as the 
reciprocal of a lower cutoff on |fc — g]. We find that {2£ + 
1) poles are generically present for each i when Fourier 
transforming 



G{k) = Y. 



2e + i 

Ait 



i 

(C4) 



which simplifies to 



knT 



-\eUk)\'. (C5) 



e - Ki;\\{k)+^,i 
Our two-dimensional soccer-ball (^ = 1) Hamiltonian: 

H = -Y^J{x)* S{y) (C6) 
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leads to three poles at all temperatures. Much unlike the 
pole locations, the branch points are temperature inde- 
pendent. If we ignore the k dependence of |ef Q(fc)p, then 
the contribution to from the branch cut contour 

will read 



e 



Attx- 



■{A+-A_) 



/ dv 

Je/A 



V ^/e^ - 167;V.x2 /o(u) 



16fi)' 
(C7) 



where A± = ^le^QpfesT and Iq a Bessel function of 
imaginary argument. The contributions from the branch 
cut would seem (if {e^^ q} were k independent) to give 
rise to quasi-algebraically damped correlations. 
For the case at hand, 



AiOky 



AiOkx 



\u+)=N+\ 



Ai9k„ 



-1), 

.1), 



are the eigenvectors expressed in the Cartesian basis. In 
the spherical tensor basis {|e±°|2} are the squared moduli 
|A^±|2, of the z-components of the eigenvectors \u±) are 



= 1,01 



= 1 - 



8k^9 



2a2 



± 6l3Vl6fc2 + 612 + IQ02},2 ■ 



(C8) 



The dependence of \e^^i on \k\ is manifest. Insert- 
ing this in G{k) and Fourier transforming one readily 

verifies that no quasi-algebraic behavior remains only 
the usual exponential correlations originating from sim- 
ple poles appear. For infinitesimal curvature (or frustra- 
tions) 0, all correlation lengths (originating from these 
poles \Im{ko}\~^ , \I'm{k±}\~^ can be made as large as 
desired when T ^ 0: kl = -[n + e% fc^ = i[6l2 - 2/u ± 

is no crossover of 

a correlation length into a modulation length: The poles 
{±fc^} are complex conjugate for /x > /Xmin = — 7^2/i5 
and will give to a damped oscillatory behavior at all tem- 
peratures. Nonetheless, a crossover does occur (albeit 
not a thermodynamic one). At very low temperatures the 
long correlations are spawned by k± and arc arc therefore 
oscillatory. At intermediate temperatures the correlation 
length generated by fco can be larger leading to asymp- 
totic long range correlations that are non-oscillatory in 
character. To see this, we note that 

\Im{ko}\ = 
|/m{fc±}| = (2^^ + Sm^^ ^ ^2)i/4| gj^(^/2)| 



, Vl6Mg2 + g4 



(C9) 



Whenever ii{T) = 0^/2, 



|7m{fco}| = 

\Im{k±}\ = (3.75)1/4 1 /TO{fco} I . ^^^q^ 

At this temperature |/m{fc±}| > |/m{fco}| and therefore 
the asymptotic correlations will not be oscillatory as they 
are near the ground state. In this way we can identify a 
low temperature crossover. It is straightforward to obtain 
a simple expression for the value of the chemical potential 
H at this crossover temperature. 

|/m{fc±}| = 2-i/2(2(?4 + 3^02 _^ ^2)1/4 



-1/2 



16/7,6*2 + 6*4 

(02 _ 2/i)2 ^ 

\Im{ko} \ = 



(Cll) 



For larger matrices (e.g. ^ > 1 in an 50(3) action) 
there are potentially more crossovers. For these matri- 
ces, at asymptotically large temperatures all correlation 
lengths iJL~^/'^. As we will see, this is not the case for 
long range interactions Notice that this change from a 
uniform behavior at high T to finite modulation lengths 
at low T is a feature shared in common with the Coulomb 
frustrated ferromagnet (and any other system with a frus- 
trating long range interaction). We find similar results 
for the 50(4) action in three dimensions of [7]. The cor- 
relator 



G{k) = Y.- 



kuT 



27r2 i^„A-=o,i(l^l)+Mn 



(C12) 



with i a single index labeling the eigenvalues, where we 
note the commutativity of [M+ — M_] with the kernel. 



2. Correlation Function Crossovers Within the 
Coulomb Frustrated Ferromagnet 

Unlike other Appendices, most of the contents of 
this section appear in a Letter version of our work on 
Coulomb frustration [11]. The continuum limit Hamilto- 
nian of the Coulomb frustrated ferromagnet 



H : 



1 /^.2 



'k-^ + ^)\S{kf 



(C13) 



where the spherical constraint (global spin normaliza- 
tion) is enforced via the chemical potential fi. A non- 
increasing weight for all modes k mandates ji > — 2e. 
The minimizing modes = e). By equipartition, the 
mode occupancy {\S{k)\'^) = ksT/iP + e'^/P + n). The 
correlator: 

G(f) = — i-^ / Sk (|S'(fc)|2) exp[zfc-f] 



(27r)3 

keT 
27r2|f| 



dk 



fc3[/m{e*'=l^l}] 

(fc2+a2)(fc2-h/^2) 



(C14) 
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where 



2 2 M T - 4e2 



When /i > 2e, the integral can be readily evaluated by 
applying the residue theorem to the poles lying on the 
imaginary axis at fc = ±ia, ±1(3, 



G{x) 



ksTiJ-c 



47r|f|(/?2 -a2) 



(C16) 



Note the existence of two macroscopic correlation lengths 
- a consequence of charge neutrality. As the spins portray 
charges, they must sum to zero, 

j G{x) d^x = {\ j S{x) d^x]"^) = 0. (C17) 

Whenever G is dominated by its long-distance behavior, 
the integral can vanish only if G{x) contains positive and 
negative contributions, as in Eqn.(C16). This integral 
can vanish only if G{x) contains, at least, two length 
scales. At high temperatures, the length 

h = \Re{P}\-'^ « for fj, » 2e). (C18) 

plays the role of the correlation length of the canonical 

short-range ferromagnet (i.e. with no frustrating chargc- 
e = 0). Note that now, however, an additional correlation 
length appears: 



h = \Re{a}\-' 



1 

1^' 



(C19) 



For /x^ = 4e^, we find that a = fi, and the four poles 
merge in pairs. When /i^ < 4e^ the poles reside on a circle 
of radius e^/"*; i?e{a^} = i?e{/3^} = |. This corresponds 
to a pole at k on the circle of radius e^/^ at an angle 
6 with cos2(? = — r/(2e). The analytic continuation for 
/i^ < 4e^ is trivially {a = ai+ ia2 =/?*): 



G{x) = ^^^^ exp{—ai\x\ 



47r 



{0.2 — Oil) sin oi2\x\ + 2a\a2 cos 0.2 
4aia2|^| 



where a = ai + ia2. 

Putting all of the pieces together: 



(C20) 



for T > Ti = fcBi(27r)V[ ^ dHji^e +^ + fc^)] 

>Te(e = 0); £ij^£2, (C21) 



Ld 



2[/x-|-2e] , and a single modulation length 



27r/|a2| = 47r[-/x + 2e] 



-1/2. 



at high temperatures 



(C15) (T > Ti), the system possesses two distinct correlation 
lengths. When T — , the modulation length diverges 



as Li) ~ (Ti — T) The crossover temperature and 

the critical temperature read 



Stt 1/2 



Tc{e = 0) = 



A ' 



(C23) 



with A an upper momentum cutoff. 

The mechanically oriented reader might note a simple 
analogy to the temporal quadratic form for the damped 
oscillator, where two relaxation times appear for 7^ > 
4mK (with 7 , K, and m the dissipation coefficient, the 
spring constant and the mass respectively). At high tem- 
peratures, ii « /i"^/^ plays the role of the usual correla- 
tion length (for a free Landau- Ginzburg action with the 
quadratic kernel u(A') = /i + fc^); this is the analogue of 
the damping time in an "un-drivcn" mechanical system. 
At low temperatures (T < Ti), the information previ- 
ously encoded in two correlation lengths at high tem- 
peratures is now manifest as a single correlation length 
e = UJ-i = 2 



Ld = 2n/\a2\ 



2^fj2 and a single modulation length: 
47r. 



2e-n 



. The "domain size" Ld, the 



inverse characteristic sinusoidal modulation length of the 

correlations, diverges at T = Ti, the correlation length at 
Tc = 0. The value registered by /j. at the avoided critical 
temperature Tc{e = 0) is roughly dimension independent 
and pinned at /i « — e; both Ld and £ are regular at 
Tc{e = 0). Note that, lime^o+Tiie) = Tc(e = 0) (this 
statement can be well defined only in a slightly modi- 
fied model such as v{k) = k^ + e^fc~^ + A ^i^j kfkj with 
small A). In that limit Ld makes an entrance at the 
"avoided critical point" Tc(e = 0). Such a crossover of 
a correlation length(s) into a modulation length(s), or 
vice- versa, is a feature that is generically absent in finite 
Range interaction systems. 



3. Correlation Functions for Other Short Range 
Frustrated Models and others 



Consider the interaction with the kernel 
v{k) = [A^ - 
where the lattice Laplacian 



(C24) 



for T < Ti £1= £2. 



(C22) 



The temperature Ti, defined by ii(T = Ti) = 2e, 
marks a dramatic crossover. At low temperatures (T < 
Ti), the system possesses a single correlation length ^ = 



A(fc) = 2^(1 -cos fci), 



(C25) 



with {ki}f_i the momentum components. As the lat- 
tice Laplacian connects, in real space, two sites that 
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are a lattice distance apart, the kernel of Eqn.(C24) 
may link sites that are, at most, six lattice constants 
apart. It corresponds to a Range = 6 interaction. If 
the constant Aq > then upon expanding Eqn.(C24), 
we will find both positive and negative contributions. 
Expanding Eqn.(C24), we find both ferromagnetic and 
anti-ferromagnetic interactions that compete on differ- 
ent ranges. When a chemical potential term is added to 
enforce normalization, we find that its minimum value is 
— ^0 (whence the mass gap vanishes). For future 
reference let us denote the mass gap by c = (/x — Umin)- 
The poles of the propagator 



G{k) = kBT[c + v{k)]- 



(C26) 



occur at 



{A±}3^i = [A3 ± zVH]i/3 exp (27rm/3). (C27) 

The poles {A^} form the vertices of a hexagon in the 
complex k plane. Whenever fe| = fc^ with ko a given 
complex number, then there must exist, at least, two 
roots with different imaginary values {\Im{ki}\). The 
correlation lengths ii = \Im{ki}\~^ . Two different cor- 
relation lengths £i ^ £2 can be identified at all /i > iJ.min 
[n = Umin only at T = Tc = 0). We note the exis- 
tence of multiple correlation lengths £i at all values of 
Ao and temperatures (the mass gap c). Multiple cor- 
relations lengths are also present everywhere- this also 
includes the unfrustrated system having Aq < 0! Mul- 
tiple correlation lengths are generically present for func- 
tions of A (or fc^ in the continuum) and are accompanied 
by a discontinuity in the critical temperature in suf- 
ficiently high dimension whenever appropriate competi- 
tion allows real roots < {Aj} < Ad when the mass gap 
vanishes (/i = jimin)- For the canonical A - polynomial, 
the system possesses several modulation lengths and sev- 
eral correlation lengths at all T; barring few exceptions 
- their net number is conserved as temperature is varied. 
Generically, all cross-over temperatures [Tj=i^2,...,p] (in- 
cluding low dimensional systems which possess no crit- 
ical behavior for zero frustration), at which correlation 
lengths disappear and turn into modulation lengths, tend 
continuously to the avoided critical temperature (or its 
analytic continuation for low dimensions - in high dimen- 
sions this temperature /x(T') = fi* becomes critical for 
zero frustration). For a general finite ranged A poly- 
nomial kernel, the system possesses a fixed number of 
correlation lengths and a fixed number of modulation 
lengths; there is no sharp analogue of crossover temper- 
atures {Ti\ wherein modulation lengths turn into corre- 
lation lengths. These crossover temperatures are more 
prevalent for non- analytic functions of A. The Yukawa 
Ferromagnet v{k) = fi+ {k'^ + X^} + -gr^ (in real space 
the last term gives rise to Yukawa like screened interac- 
tions) has Tc(e = A^) > in d > 4 and in any dimen- 
sion Tcie > A^) = 0. Two correlation lengths appear 



for /i^ > 4e^ (including all unfrustrated ferromagnets 
(e<0)). 

For the non-analytic kernel v{k) = {\k\ — q)^, 



-tan ^(g/V^)} + 9ln( 



M + (A - qf 
H + q^ 



)], (C28) 



with Umin =0. 

The correlator corresponding to v{k) = {\k\ — q)"^, 



G{x) = ^ ^.^. e ^1^1 {q sin qx -\- y/Ji cos qx] 
2^\x\ 

ql{\x\ 



(C29) 



where 



f 

Jo 



dn 



«;exp[— K|a;| 



{iJi + q^ — K^Y + Aq'^K^ 



> 0. (C30) 



Note the algebraic decay in this artificial non-analytic 
system. Here the system displays avoided critical behav- 
ior and exhibits one (two) modulation lengths 211/ q (and 
infinity), which does not vary with temperature. Within 
the spherical model, this behavior is an exception: mod- 
ulation lengths tend to vary with temperature. 

With t = {ksTciq = 0))~^ - (/cbT)"!), we find at 
low temperatures (/x^/^ terms neglected) for the v(k) = 
{\k\ — q)^ action. 



i« (1 + 



1 



1 



2^A^q^^ + i^ 

L 1 . M+(A-g)^ . 1 q^ 

2tt^'^'^^ 11 + q' ' 2tt^' 



(C31) 



In the low temperature /z — > 0+ limit: f ~ 

in any dimension d ^ 3 apart from a constant multiplica- 
tive factor). The large n critical exponent Vq^o = 1. At 
very low temperatures {t < with large |t|), 

V^~aVN~a' [Vo(-i)]- (C32) 

If, by analogy to the standard quadratic part of Landau- 
Ginzburg actions, we define a length Ld ~ r^^^'^, the 
inverse square root of the mass gap which is the cor- 
relation length in the canonical (q = 0) case, then as 
T 0+,Ld ~ e-if^io (with q = e^/^). We point here 
to the strong resemblance to the previously obtained re- 
sult for the exact correlation lengths appearing for the 
Coulomb frustrated ferromagnet with v{k) = [e'^/k-^ + k^] 
at high temperatures. 

For the short-range (Teubner-Strey) correlator 



G-\k) = a2 + cik'^ + C2k^ 



(033) 
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it is a simple matter to show that 

sin KX 



G{x) 



KX 



■exp[-a;/^], 



(C34) 



where 




(C35) 



In amphiphihc systems 02 and ci,2 are functions of 
amphiphile concentration (as well as temperature). In 
many simple thcrmodynamical models, short range in- 
teractions interactions thwarted by others, the modu- 
lation length always varies with temperature: in systems 
with long-range interactions, the modulation length in- 
creases as temperature is raised while the converse gener- 
ically holds for short-range interactions. 

When the kernel v{k) = k^, the real space correlator: 



G(x) = 



47rXy//I 



exp [— a;/^] sin kx. 



(C36) 



Within the spherical model, this ferromagnetic system 
displays thermally induced oscillations. At T = the 
(ferromagnetic) ground state is unmodulated. 



APPENDIX D: THERMODYNAMIC FUNCTIONS 
IN THE VICINITY OF THE DISORDER LINE 

1. Non-Abelian Backgrounds 

For a global gauge field normalization, all multipli- 
ers {/x^}"^]^ = const = II,. The internal energy density 
U /N = n{kBT — n). AS before, the chemical potential n 
may obtained by inverting l/lksT] = J -^^G{k). For 
fixed (non-dynamic) non-Abelian backgrounds, there are 
no explicit high temperature crossovers- we find no ther- 
modynamic signatures of a disorder line. The number 
of correlation lengths and modulation lengths is temper- 
ature independent. However, the long range physics is 
sinusoidal (governed by the poles k±) at low tempera- 
tures while it is nonsinuiodal (dominated by the pole at 
ko) at intermediate temperature. 

We now examine more complicated gauge theories of 
glasses. To low orders in perturbation theory, our con- 
clusions are not different. In the ^— th representation of 
50(3) model, the eigenvalues for k6 ^ 1 are 



Xm{\k\)^e^[ik9f +2m,ik9) 
+l{l{l + l)+m'} 

1 

[{I + m){l - m + 1) 



x{l + m- l){l-m + 2) 
-{I - m){l -I- m + 1) 
x{l + m + 2){l-m + l)\]. 



To this order, the equations 

M + A™(|fc|) = 



(Dl) 



(D2) 



have a pair of complex conjugate zeros and a single real 

zero when /i > /imin- There is no crossover of a pure imag- 
inary (fc) poles into complex ones. The computations are 
similar for an 50(4) background gauge. When kO <^ 1, 
the eigenvalues may be expanded as polynomials in (kO). 
The behavior of the correlation lengths and thermody- 
namic functions is the same as that observed in the frus- 
trated short range systems that we have discussed no 
merger of poles occurs, the number of correlation lengths 
and the number of modulation lengths are temperature 
independent. Thus fixed non-Abelian backgrounds do 
not lead to explicit analytical (continuation) crossover in 
the thermodynamical functions. 



2. Thermodynamics near the Disorder Line for The 
Coulomb Frustrated Ferromagnet 

In the spherical (n — > 00) limit, the internal energy 



_ J_ \ - fee 



OM^ ^^v{k)^\{kBT-„). (D3) 
2A ^ v{k) +11 2 



The chemical potential, /z, may be computed as a func- 
tion of the temperature T by inverting the spherical equa- 
tion of state -rr-Tr = — S°(r = u) with YP the zeroth or- 
der self-energy. YP{r = /x) attains different forms for 
> 4e2 (when T > Ti) and < ie^ (T < Ti). 
The relation tan~^ ^ ~ ^ applied to re- 

late the two natural functional forms. Although there 
are no divergences at the merger point of the poles 
T = Ti, the explicit functional forms of the suscepti- 
bility ilSigjl"^} = and the specific heat Cv = ^ 

do change. For T > Tc, ^{T) is monotonically increasing 
in T. 



1 



rmip = 0) ' 



(D4) 



with the zero momentum polarization diagram n(p 



0) = 



J {v{k)+ixy 



Here we see that no specific 



32[2{ke) + m + 1] 



heat extrema can be encountered at finite temperatures 
in any spherical model for fi > /imin {T > Tc). In 
limn 00 (and as shown, also to 0(l/n)), Ti is an an- 
alytic crossover temperature. As many thermodynamic 
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entities are integrals of pair correlators, the change in 
G{x) above and below Ti suggests a crossover in cer- 
tain thermodynamic quantities. We expect this crossover 
to persist to all orders in 1/n. We indeed find this to 
be the case to 0(l/n). The underlying reason for this 
crossover is very transparent: In the factored propagators 
+ oe^){k^ + tbe poles have different explicit 
forms, if T is above or below the disorder line temper- 
ature Ti. The critical temperature Tc{e), as computed 
within the spherical model, drops discontinuously. Above 
Ti we have two correlation lengths. Below Ti there ex- 
ists a single modulation length and a single correlation 
length. As we have just outlined, even though true long 
range order does not set it until T < Tc(e), an analytical 
crossover occurs at T = Ti(e). At low temperatures this 
crossover temperature becomes the avoided critical point 
Tc(0). We have suggested that Ti can be regarded as 
an effective dimensional crossover temperature. At lower 
temperatiires the behavior becomes more and more like 
that of a one dimensional system. At higher tempera- 
tures, the behavior is, in most respects, similar to that 
of a three dimensional nearest neighbor ferromagnet. 



APPENDIX E: A GENERALIZED 
MERMIN- WAGNER INEQUALITY 

1. The Classical Case 

Our approach is the standard one. We will keep it more 
general instead of specializing to anti/fcrromagnctic or- 
der or to interactions of one special sort. With the nota- 
tions introduced in Section(II), we investigate n compo- 
nent spins on a lattice. An applied magnetic field 



h{x) = hcos{q ■ x)ec 



(El) 



causes the spins to take on their ground state values. 

If n = a = 2 the unique spiral ground state {S^{x)), to 
which a low temperature system would collapse to under 
the influence of such a perturbation is 

Sf{x) = sin{q- x). (E2) 

When n = 3 the ground state is not unique: 

Si<ni^) = '^^ sin(g-f) 

n-l 



and a magnetic field may be applied along two directions, 
with all the ensuing steps trivially modified. With the 
magnetic field applied 

n 

^ = 2 5 ^ ~ - E Kix)Sn{x). (E4) 

x,y i=l X 



Note that the knowledge of the ground state is not imper- 
ative in providing the forthcoming proof [30] . We exploit 
the standard rotational invariance of the measure 

J dn ■= Z-^ J llrS{x)6(S\x) - l)e-^". (E5) 

X 

Note that this is not applicable to many other spin mod- 
els (e.g. Dzyaloshinskii-Moriya interactions). The gener- 
ators of rotation in the [af3] plane are 



-a/3 - '^aggff 

d 



d 



(E6) 



= -Jd-S6iS^-l) 

f{Si, Sa cos6 + Sj3 sinO, ... 
, S/} cos6 — Sa sin a, Sn). 

= j d"S5{S^ - l)Laf3f{S). 



(E7) 



(E8) 



In the up and coming, _L will denote the the projection 
along the /? direction. We define the operators 

A{k) = y~^exp[ifc • x]S±{x), 

X 

B(k) = Eexp[i(fc g) • x\Ls{(3H). (E9) 

X 

By the Schwarz inequality 
K E 4*A>P < ( E A*^.) * ( E B*B^). (ElO) 

We will let i = a, /? in the sum span a two element subset 
of the n spin components. For any functional C: 

Lsie-^"C) = e-^"{Ls{C) + CLs{-(3H)}- (EH) 
° = / n '^"^(^) <5(^2(x) - 1) Ls[e-l^"C] (E12) 



{CB{p)) = {Y,^Mip-ALx{C)) 



(E13) 



E(4(4(/3^))) =/3([E Si{x)Si{y)V{x-y) 

i=a^j3 i=o:,j3 

-h{x)Sn{x)]) (E14) 
{B{kr ■ B{k)) = /3E{(cos(^ + q)-{x-y)-l) 
[Y,{SmSM)\v{x-y)} 

i=a,l3 

-hix){Sn{x)} > (E15) 
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Henceforth, for simplicity, we specialize to n = 2. 
Fourier expanding the interaction kernel 



V{x-y) = -J2''{r)e'''^^-y\ 



(E16) 



and substituting 

{S{x) ■ S{y)) = ^(|5(n)|2)e™-(-^-«^, (E17) 



we obtain 



< {B{ky ■ B{k)) = 0AfE - ph^{Sni-q)) = 

/3 



n 

-2v{u)] - (ih^{Sn{-q)) (E18) 



(2) 

where A-^ E measures the finite difference of the internal 

k 

energy with respect to a boost of momentum k. 

{A{kr-Aik)} = Y,{s±{x)-s^m 

x,y 

X exp[ik ■ {x - y)]. (E19) 
{A{kr ■ B{k)) = (^Lt{S^{x)) exp[j(fc + g) • x\) = m^, 

i,x 

where = {Sn{q)) and, as noted earUer, _L refers to 
the i = 1 spin direction orthogonal to « = n = 2. Note 
that with our convention for the Fourier transformations, 
a macroscopically modulated state of wave-vector q, the 
magnetization niq = 0{N) as is the energy difference 
in Eqn.(E18). Trivially rewriting the Schwarz inequality 
and summing over all momenta k, 



^ {Bikr . B{k)) - Y 

which explicitly reads 

2N\m,^\^[0^i{\Siu)\'Mk + u) 



(E20) 



+v{u ~ k) - 2v{u)] + 2\ 



r 



k ^'V 

= Nj20lm- (E21) 

X 

Explicitly, as the integral j''"''"^ non-negative 
(as (Bik)* ■ B{k)) > the denominator in Eqn.(E21) is 



positive for each individual value of fc), and as {S']_{x)) < 
1, we obtain in the thermodynamic limit 



\k\<S 



{2ttY 



{\S{u)\'^){v{k + u) + v{k - u) - 2w(m)) 



+2\h\\mq\ 



< 1. 



(E22) 



Taking S to be small we may bound from above (for each 
value of k) the positive denominator in the square brack- 
ets and consequently 



r\k\<S jd^. r ^ 
J (27r)4y {2w\ 



{\S{u)n{v{k + u) + v{k-u) 



+2\h\\7 



1 ^ r\k\<5 ^d^ r dSl^ 
-J Wi^VJ (2^ 



u 

Y 

2v{u)) 



A,k^\s{\Siu)\^) + 2\h\\mg 



(E23) 



with A^j chosen to be the largest principal eigenvalue of 

the d X d matrix didj[v(u)], and Ai a constant. For a 
twice differentiable v{u), and for |fc| < ^ where S is finite, 

{v{k + u)+ v{k -u)- 2v{u)) < AiX^k'^ < Bik'^ (E24) 

for all u within the Brillouin Zone with the additional 
positive constant Bi introduced [31]. Here we reiterate 
that {B*{k) ■ B{k)) of Eqn.(E18) is positive definite and 
consequently the bound derived is powerful. 

In d < 2, the integral j''^'^'^ diverges making 

it possible to satisfy eqn.(E22), at finite temperatures, 
when the external magnetic field /i — > only if the mag- 
netization niq = 0. If finite size effects arc restored, in 
a system of size N = L x L where the infrared cutoff 
in the integral is 0{^,^) the latter integral diverges 
as C'(lniV). This implies that the upper bound on |m^| 
scales as 0{N/^lnN), much lower than the 0{N) req- 
uisite for finite on-site magnetization. For further details 
sec [30]. If there are M > 2 pairs of minimizing modes 
and 2p + 1 > n > 2p (with an integer p) then we may 
apply an infinitesimal symmetry breaking magnetic field 
along, at most, min{p, M} independent spin directions 
(a). Employing the spin rotational invariance within 
each plane [a/3] associated with any individual mode £ 
(au lieu of a specific q) wc; may produce a bound similar 
that in Eqn.(E22) wherein (|5'('u)|^) will be replaced by 
Ei=a,/3(l^i(^)H and |m,H ^ |m„(/)|. 



2. The Quantum Case 

The finite temperature behavior of a quantum system 
is, in many respects, similar to that of a classical system. 
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The quantum system is also invariant under rotations 
with S^{x)=S{S+l). 

Alternatively, one could directly tackle the n = 3 quan- 
tum case by applying the Bogoliubov inequality 



{{A,A^).{[[C,H],C^])>\{[C,A])\^ (E25) 



with [ , ] and { , } the commutator and anticommu- 
tator respectively. Prom this inequality, it follows that 
([ [C, H],C^) is positive definite. In particular, for any 
six operators {^i, A2, A3, Ci, C2, C3}, 



0=1 



J2m,Ai])f (E26) 



Setting Ai = S2{q — k) and Ci = Si{k) we will once 
again obtain Eqn(E22) with the classical spins replaced 
by their quantum counterparts. 

Rather explicitly, employing 

[S°'{k),S^{k')] = ie°'^"'S^{k + k'), (E27) 
we find for the Hamiltonian of Eqn.(E4) (with n = 3), 

[[Ci,H],Cl] = -LY^^S2{k')S2{-k') + Ss{k')Ss{-k')) 



X [v{k + k') - 2v{k) + v{k' - k)] 
+ lT.hs{k')Ss{-k'). 



Similarly, 



\{[C^,A^])f = \{S,iq))f, 



(E28) 



the squared magnetization along the z (or 3) direction 
for a mode q, and 



X^({Ai,At})=2^(52(fc)52(-fc)) 
fe k 

= 2{[S2{S=0)f). 



(E29) 



Next, let us cyclically set, A2 = 5*3(5 — k),C2 = 
S2{k),A3 = Si{q-k), and C3 = Ssik). The commu- 
tators associated with these operators are all identically 
the same apart from a uniform cyclic permutation of all 
spin components involved. 

Trivially rewriting the symmetrized Bogoliubov in- 
equality Eqn.(E26) and summing over all modes k, 



k,a 



(E30) 



Replacing the k sums by integrals in the ther- 
modynamic limit, and employing the positivity of 
([[Cq, H],Cl]) that follows from the Bogoliubov inequal- 
ity for each individual value of k, we find the trivial quan- 
tum analogue of Eqn.(E22), 



1 

2P 



\k\<5 



{2n) 



{\S{u)\^){v{k + u)+ v{k -u)- 2v{v)) 



+ 



m„ 



<S{S+l). 



(E31) 



Apart from the simple scaling factor of AS{S + 1) = 
2 X]a({Ao, A\}) by comparison to the classical case, there 
is no difference between this inequality and its classical 
counterpart in the zero field limit. We symmetrized the 
Bogoliubov inequality (Eqn.(E26)) in order to avoid the 
appearance of only two transverse spin components in 
the Fourier weights {\Si{q)\'^) so as to give the resulting 
Mermin- Wagner inequality a transparent physical mean- 
ing associated with the energy boost differences ' E 
which are symmetric in all spin indices. From here on- 
ward the discussion can proceed as in the classical case. 



APPENDIX F: MERMIN- WAGNER BOUNDS IN 
HIGH DIMENSIONS 

In any dimension, Eqn.(E22) reads in the limit /i ^ 



2\mffVT 



(fik 



(2^Y ^(2)^ 
k 



< 1 



(Fl) 



with the shorthand defined by Eqn.(E18). By parity in- 
variance and noting that {S{x)} are real, 

Et^(fc + «)(l^(u)|^) =E^(^-^)(I^(^)I')- (F2) 



Thus the denominator of Eqn.(Fl) reads IS.\'E = 
[2(£'g — Eq)] where E-^ is the internal energy of system 
after undergoing a boost of momentum k and Eq denotes 
the internal energy of the un-boosted system. In some 
instances when the dispersion relation about an assumed 
zero temperature ground state is inserted into Eqn.(Fl), 
we will find that the integral in Eqn.(Fl) diverges: At 
arbitrarily low temperatures we cannot assume the zero 
temperature ground state with the natural dispersion re- 
lation Ai? for fluctuations about it. A case in point is the 
dispersion relation for the Coulomb Frustrated Ferromag- 
net. The denominator in Eqn.(E22) is a finite tempera- 
ture extension of the T = dispersion relations (A^ 'E). 
In general, at zero temperature. 



{\S{kt) = ^[6j:^^- 



(F3) 
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and the integral in Eqn.(Fl) becomes Eqn.(35). 

Whenever Eqn.(35) diverges, an assumption of an al- 
most ordered ground state at arbitrarily low tempera- 
tures (T = 0^) with the ensuing zero temperature disper- 
sion relation {Al ' E — [^{E^ — Eq)]) about that ordered 
state is flawed: Eqn.(Fl) is strongly violated. This poses 
no problem for most canonical d > 2 dimensional models 
where the integral in Eqn.(35) is finite. For the three 
dimensional Coulomb Frustrated Ferromagnet and other 
high dimensional models the divergence of this integral 
hints possible non-trivialities. The average decoherence 
(or relaxation) time scale (r) may diverge if we assume 
a zero temperature dispersion of fluctuations about an 
ideal crystal. A divergent decoherence time scale suggests 
the absence of broken translational symmetry. Formally, 



/ 



(27r)3 



Tk = T{f= 0). 



(F4) 



In order for magnetization = 0{N) to arise, the av- 
erage of the inverse boost energy over all of k space 

1 



>r < 0{T-'), 



(F5) 



decoherence) time of the system averaged over magnons 
of all possible momenta. Whenever the average charac- 
teristic relaxation time 



r d^k r r 
J (27r)4y 



d'^i 



{\S{u)\^){v{k + u) + v{k-u) 



(27r)<« 
-1 -1 

2v{u)) 



diverges then by our generalized inequality 
2|m<,|2T(rg) = 2\mg\^TTr=o < 1 



(F7) 



(F8) 



the system does not order in such a way that the fluctu- 
ation dispersion about any viable ground state is valid. 
Moreover, as we show in the text the system displays 
divergent glassy dynamics. The characteristic divergent 
relaxation times in glassy systems suggest a divergent de- 
coherence time and an inability to write a wave function 
of a quantum glass. The glassy state may be specified by 
a density matrix [36]. 

* Present address: Theoretical Division, Los Alamos 
National Laboratory, Los Alamos, NM 87545 



at all < r < Tc. In most d> 2 systems this is trivially 
satisfied with the average bounded by a constant at zero 
temperature. This average diverges at T = whenever 
Eqn.(35) does. There are two possibilities: 

(i) The system is ordered at all temperatures T < 

in which case, the thermodynamic average of Eqn.(F5) is 
finite for all T > and r is non-analytic at T = 0. 

(ii) The system is disordered at all finite temperatures 
and orders classically only at T = 0. 

The first possibility ((i)) was argued for by a non- 
rigorous yet elegant diagrammatic analysis by Brazovskii 
[35] long ago: Thermal fluctuations, on their own, may 
enhance (or generate) cubic terms fortifying (or trigger- 
ing weak) first order transitions. This cannot be ruled out 
by the rigorous Mermin- Wagner inequalities that we de- 
rived. As reiternated, if the fluctuation integral diverges 
then we may not obtain the low temperature T = 0+ 
dispersion by assuming a nearly perfectly ordered state. 
This does not rigorously preclude order. Order, if it exists 
at arbitrarily low temperatures must display non-trivial 
excitation spectra about it. If such a possibility arises, it 
might be immaterial if the system is permanently frozen 
into a glass before reaching an equilibrium thermody- 
namic transition. In both cases, t{T) is non-analytic 
at T = 0. The integral of Eqs.(35,F5) has a suggestive 
physical interpretation. If the quantum spin system is 
subjected to a boost of momentum k, then 



1/AfE ^ r, 



(F6) 



is the characteristic lifetime of the excited state. The 
average in Eqn.(F5) is the characteristic relaxation (or 
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